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ABSTRACT 

. Magnetic elements, inductors and/or transformers are among the most 
important elements found in electronic power processing and filter circuits. 

In filter circuits, the inductor is used to smooth the current flowing to 

i 

the load and its output voltage. In energy-storage dc-to-dc converters 
used in power processing applications, the inductor or transformer is used 

1 

to store and transfer energy in some controlled manner between source and ' 

load. The magnetic element is nearly always a custom designed and con- . 

i 

I 

structed element as opposed to other elements used in these circuits. There- ! 
fore, the design of the magnetic element is a crucial step in the design of l 

power processing and filter circuits as a whole. In high-power applica- [ 

tions, usually it is desirable to design these inductors using air-gapped 
magnetic-core structures. This dissertation presents methodical approaches 

I 

to the design of inductors for use in LC-filters and d‘c-to-dc converters 
I using air-gapped magnetic structures. 

I In addition to the design procedures for the inductors, methods for 

the analysis and design of full -wave rectifier LC-filter circuits operating 
with the inductor current in both the continuous-conduction and the discon- 
tinuous-conduction modes are presented. In the continuous-conduction .mode, ' 

I 

I 

(iii) 

I 
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linear circuit analysis techniques are employed, while in the case of the 
discontinuous mode, the method of analysis requires computer solutions of 
the pi ecewise-1 inear differential equations which describe the filter in 
the time domain. Using these methods of analysis, a family of design curves 
is obtained relating filter parameter values of inductance, capacitance, and 
load resistance to the performance factors of percentage output ripple 
voltage, average output voltage level, rms and peak transformer-winding 
currents, and input power factor. 

Three procedures for designing filter inductors using air-gapped cores 
are presented. The first procedure requires digital computation to yield a 
design which is optimized in the sense of minimum core volume and minimum 
number of turns. The second procedure does not yield an optimized design 
as defined above, but the design can be done by hand calculations or with 
a small calculator. The third procedure is based on the use of specially 
prepared magnetic core data and provides an easy way to quickly reach a 
workable design. 

The procedure for designing air-gapped energy-storage reactors for 
dc-to-dc converters is based on an earlier study of the storage and transfer 
of energy in the magnetic reactors. The continuation of a study in this di- 
rection leads to a simple relationship for the required minimum volume of 
the air gap. Determination of this minimum air-gap volume then permits the 
selection of either an air-gap length or a cross-sectional core area, thus 

identifying a workable magnetic structure. 

It is believed that the analytical guidance in conjunction with the 
design procedures presented in this dissertation can significantly reduce 


(Tv) 



the time, and hence the cost, of designing full -wave rectifier LC-filter cir- 
cuits, air-gapped magnetic-core inductors for use in filters, and energy- 
storage reactors in dc-to-dc converters using air-gapped magnetic structures. 
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AIR-GAPPED STRUCTURES AS MAGNETIC ELEMENTS 
FOR USE IN POWER PROCESSING SYSTEMS 


j 


CHAPTER I 
INTRODUCTION 

One of the most important elements found in electronic povjer processing 
and filter circuits is the magnetfc element which may be an inductor and/or a 
transformer. While other components in these circuits, such as semiconductor 
elements and capacitors, are chosen from populations of standard sizes and 
ratings, the magnetic element is nearly always a custom designed and constructed 
element. As power processing circuits and systems have become more complex, the 
design of the magnetic element has necessarily become a more difficult and subtle 
art. Few elements available to the power processing engineer can be designed 
so closely to the desired parameters as the magnetic element, even though there 
are standard size and permeability groupings available for the magnetic cores 
on which windings are placed. While giving the designer a greater range of 
choice in circuit planning, this extra degree of freedom is not without its 
disadvantages. The complexity of the designer's problem is only partially in- 
dicated by the multiplicity of nonlinearities in the large range of magnetic 

materials available. The variety of structures and shapes is extensive as 
well, with shapes ranging from presized toroidal cores with known effective 
permeabilities to cut and pot cores for which the designer must determine the 


( 2 ) 
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air gap and, therefore the effective permeability. In every case, the primary 
goal of the designer is to determine the proper number of turns and wire size 
required, given that a particular magnetic structure has been selected, so that 
some set of design parameters will be achieved. While the design parameters 
frequently are a specified inductance value for a given dc current, there are 
other important criteria in power processing circuits which lead to different 
specified design parameters. By approaching the inductor design problem from 
a broader point of view, it is possible to take into account other factors of 
great importance in overall system design, such as size and weight. 

In an inductor-capacitor (LC) filter circuit which follows a full-wave 
rectifier, the inductor is used to smooth the current flowing to the load and 
capacitor due to the rectifier and its output voltage v^hich consists of uni- 
directional half sinusoids, thus producing across the load a dc output voltage 
with a small ac ripple superimposed. The filter can operate in one of two 
steady-state modes, usually referred to as the continuous-conduction and the 
discontinuous-conduction modes. In the continuous mode, current through the 
filter inductor L is always greater than zero. In the discontinuous mode, 
current through the inductor is equal to zero during a portion of each steady- 
state cycle. Various procedures have been reported in the literature for the 
analysis of rectifier-filter circuits; hov/ever, most of the design techniques 
employed today are based on procedures developed in the 1930's and 40 's and are 
limited to considering RC filters, or LC filters restricted to the continuous- 
conduction mode [1-5]. Recently, an analysis of a rectifier with LC-filter 
circuit operating in the discontinuous mode was presented under the assumption 
that the value of the capacitance C was so large that the output ripple voltage 
was negligibly small [6]. 



4 


In energy-storage dc-to-dc converter circuits, the inductor serves as 
a medium of exchange or transfer of energy between source and load. In the 
past, design procedures for the energy-storage reactor have been developed 
particularly with powdered permalloy and ferrite toroidal cores in mind [7-12]. 
These procedures make use of a collection of data for a finite population of 
core sizes and discrete effective permeabilities. In high-power applications, 
it may be advantageous to design and fabricate these inductors using air-gapped 
magnetic cores. The use of air-gapped cores permits the selection to be made 
from an innumerable population of core structure-size-permeability combinations 
and makes available a wider range of magnetic core sizes, particularly in the 
larger volume sizes. 

An important characteristic of inductors for these applications is the 
ability to carry direct current without significant change in inductance over 
a specified range of circuit operating conditions. To avoid saturation of the 
magnetic material due to the dc component of winding current, it is necessary 
to introduce an air gap in the magnetic material. In the powder cores, the 
powdered magnetic material is insulated with a ceramic material which effectively 
provides a uniformly distributed air gap. In high power applications,* it may 
be desirable to design these inductors using air-gapped magnetic-core structures. 
However, there are some disadvantages associated v/ith the use of air-gapped 
structures. Introduction of an air gap in the magnetic circuit produces couplinc 
to the fields due to other components as well as its own external fields. The 
air-gapped structures also tend to be more noisy acoustically than those without 
an air gap. 

The purpose of this dissertation is three-fold; 

(1) to present the analysis of an LC-filter driven by a full-wave rectifier 
circuit v;hich leads to a design procedure for the inductance value for the 
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j general case in which the effective load resistance Rj^ and filter capaci- 

tance C are constraint parameters and in which both discontinuous- and 
continuous-conduction modes are considered. 

(2) to develop design procedures for air-gapped magnetic-core inductors for 
use in LC-fi Iters for full -wave rectifier supplies. 

(3) to develop procedures for designing the energy-storage reactor for 
dc-to-dc converters using air-gapped magnetic structures. 

This dissertation is based in part on material presented in a series of publi- 
cations [13, 14, and 15] of which the writer is the principal author. 

I 

I Chapter II presents the analysis of LC-filter circuits with zero source 

impedance which leads to a set of curves and approximate relationships which 
are useful for designing these filters. A similar approach is followed in 
Chapter III to generate design curves for filters with a non-zero source im- 
pedance. This impedance is always present in an actual physical system, and 
its inclusion in the analysis and in the design algorithm yields a more accu- 

1 

I rate and complete design. 

j In Chapter IV, two procedures for designing filter inductors using air- 

! 

i gapped cores are presented. The first procedure requires a digital computation 
I to yield a design which is optimized in the sense of minimum core volume and 
minimum number of turns. The second procedure does not yield an optimized de- 
sign as defined above, but the necessary computations may readily be carried 
I out on an electronic pocket calculator. Using the results of the analysis 
developed in Chapter IV, and following the procedure presented by Hanna [16], 
special magnetic core curves are generated. Using these design curves, a new 
procedure for designing air-gapped inductors is presented in Chapter V. 

In Chapter VI, a procedure for designing energy-storage reactors for 
! dc-to-dc converters using air-gapped cores is presented. This procedure is 
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based on the design relationships developed in References [9—12] for gapless 
powdered permalloy and ferrite toroidal cores. 

Conclusions which can be drawn from this work are presented in Chapter 
VII with suggestions for future research. The derivation of some of the more 
important analytical expressions presented in the text are outlined in the 
Appendices. The International System of Units (SI) is used throughout, and a 
glossary of symbols and their associated units is placed at the end of this 
dissertation. 



CHAPTER II 


ANALYSIS AND DESIGN OF FULL-WAVE RECTIFIER LC FILTERS 
WITH CONTINUOUS OR DISCONTINUOUS INDUCTOR CURRENT 

2.1 Introduction 

Conversion of ac sinusoidal voltage to dc voltage for power applications 
by rectification often is accomplished using the familiar center-tapped trans- 
former rectifier or the full-wave bridge rectifier followed by an inductor- 
input low-pass LC filter which supplies the load R^ as shown in Figure 2.1. 

The sinusoidal source voltage is converted to a train of unidirectional half- 
sinusoids by the rectifier and is then smoothed by the LC-filter, producing 
across the load Rj^ a dc output voltage with a small ac ripple superimposed. 

The design of the filter for a full -wave rectifier power supply involves the 
determination of values for the inductor L and the capacitor C, the selection 
of suitable diodes, and the sizing of the transformer windings to meet the 
various application requirements. The requirements often include accomodation 
of the network to conditions of varying input source voltage and output load. 
Possible design specifications include the dc output voltage level, maximum 
allowable peak-to-peak output ripple voltage as a percentage of the dc output 
voltage, maximum allowable peak or rms currents in the diodes and transformer 
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windings, and minimum permissible power factor at the input terminals to the 
network . 

This chapter presents an analysis and a design procedure for the 
LC-filter circuit operating either in the continuous and/or in the discontinuous 
mode. In this chapter it is assumed in the analysis that the impedance of the 
sinusoidal source is negligibly small and is not included in the development - 
of design and performance relationships. However, an analysis of the LC- 
filter circuit with non-zero sinusoidal source impedance is presented in the ' 
next chapter. 

A critical condition is established as the boundary between the two 

modes of operation. For each value of the dimensionless angular frequency param 

eter defined as to/tOg, where to is the angular frequency of the sinusoidal 

input voltage source in radians per second and tOg = 1/^ is the undamped 

natural frequency of the low-pass filter, there exists another dimensionless 

parameter k, defined as coL/R^^, for which a critical value exists. For 

K > sc , the filter operates in the continuous-conduction mode; for k < k , 

Cl c r* 

it operates in the discontinuous mode. For values of k > k , the filter 

cr 

circuit is analyzed using Fourier series and linear circuit analysis techniques. 
For values of k < the filter circuit, described analytically by piecevdse- 
1 inear differential equations, is analyzed with the help of the digital compu- 
ter. An algorithm for rapid determination of the periodic steady-state condi- 
tions is used to reduce the computation time for calculations involving the 
discontinuous mode of operation. From the analyses, a 'set of curves useful for 
design purposes is generated. The derivations of useful expressions are pre- 
sented in Appendices A and B. 
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2.2 Circuit Analysis 

In this section, the rectifier-filter circuit is analyzed to obtain re- 
lationships useful for design purposes. In the continuous-conduction mode, 
linear circuit analysis techniques are employed, while in the case of the 
discontinuous mode, the method of analysis requires computer solutions of 
the pi ecewise-1 inear differential equations. In both cases, all of the elements 
of the circuit are assumed to be ideal with the equivalent load as the only 
dissipative element, 

2.2,1 Continuous -Conduction Mode 

As mentioned earlier, in the continuous-conduction mode, current ij^ 
in the filter inductor L is always greater than zero so that the source of 
voltage always is connected through forward conducting diodes to the filter 
input terminals. Assuming ideal diodes, the equivalent circuit for both of the 
configurations in Figure 2.1 is given in Figure 2.2 along with typical waveforms 
of inductor current i^ and output voltage Vq. Although the effect of diode 
forward voltage drop can be important in rectifier-filter circuits, particu- 
larly in those providing low output voltages, it is not included in the analysis 
since this factor would add an additional parameter to the family of curves 
presented later for design and performance evaluation. In one of the examples 
presented in a later section however, it is shown how diode voltage drops may 
be taken into account. The effectual value of the voltage source in the equiva- 
lent circuit, designated as vj„ is the full -wave rectified value of the sinusoi- 
dal input voltage Vj. This linear equivalent circuit can be analyzed for 
steady-state solutions using Fourier-series techniques. The effectual source 
voltage Vj is represented by its Fourier-series expansion consisting of a dc 





Figure 2.2. (A) Equivalent Diagram for the Rectifier-Filter 

Circuits Operating in the Continuous- 
Conduction Mode. 


(B) Associated Voltage and Current Waveforms from 
Computer Solutions of Circuit Equations. 
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component followed by a series of ac sinusoidal terms with amplitudes and 
frequencies given by the expression: 

wi _ /I cos2e cos4e cos6e cos2ne \ /o 

Vj “ Vj^isinei - T5 

4n -1 

where e = wt = 2-irft, where f is the ac supply frequency, and n is a positive 
integer. 

All current and voltage components in the circuit -are computed separate- 
ly by sinusoidal steady-state methods. The total current or voltage in any 
branch of the circuit may be found by summing the separate harmonic components 
from the individual terms of (2.1) using the superposition principle. Analysis 
of the filter circuit shows that its transfer function falls off at the rate 
of 40 dB .per decade (12 dB per octave) as an asymptotic slope limit. In addi- 
tion to the increasing attenuation of the filter with increasing frequency, the 
magnitudes of the sinusoidal components of the input voltage Vj also fall off 
rapidly with harmonic number. Assuming that all harmonic terms, beginning 
with the second term in the expression for Vj given by (2.1), lie in the range 
of the 40 dB per decade slope of the filter transfer function, the combination 
of attenuation and harmonic-term magnitudes yield values of 0.05 and 0.0053 
for the ratios of the fourth-harmonic and sixth-harmonic terms of Vg to the 
second harmonic term of the output voltage. Thus, the second harmonic of the 
input ac supply is seen to be the only component of significance in determining 
the output ripple voltage when the circuit is operated, in the continuous-conduc- 
tion mode, and the equivalent input voltage Vj may be approximated by the ex- 
pression: 

,.i _ ^^Im /I cos2e ^ 

IT ^2 " 3 ^ 


( 2 . 2 ) 
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The transfer function of the filter circuit operating in the continuous mode 
is given by the follovn'ng 


1 


[s^LC + s(L/Rj_) + 1] 


(2.3) 


Using the expression in (2.2) and the transfer function of the filter 
given by (2.3), a number of useful relations for various voltage and current 
components may be derived: 


DC output voltage: Vq 
P eak-to-peak ripple voltage: V^p 
Rms inductor current: I„ 


Peak inductor current: 



ir 


(2.4) 




2 4 2 

8uj^+l 6(D|^+4 k: 


1 1+(4uJ/k:^) 

^ 18(1 “8w^+16a)^+4!c^) 

1+(4cu^/ic^) 



2 4 2 

1 -8(D|^+1 6o)j^+4k 


(2.5) 

( 2 . 6 ) 

(2.7) 


2.2.2 Discontinuous-Conduction Mode 

In the discontinuous mode the action of the diodes is such that the 
current i^ is equal to zero during a portion of each steady-state cycle. The 
equivalent circuit when diodes are conducting is the same as for the continuous 
mode and is shown in Figure 2.3(A) where the second subscript 1 is added to 
signify the voltages and currents associated with this 'portion of the cycle. 

The output voltage Vq-j is obtained from solution of the following equation: . 






(C) 


Figure 2.3. Equivalent Diagrams for the Rectifier-Filter Circuits 
Operating in the Discontinuous-Conduction Mode During 
Interval (A) When Diodes are Conducting and (B) When 
Diodes are not Conducting. (C) Associated Voltage 
and Current Waveforms from Computer Solutions of 
Circuit Equations. 
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2 '* ''Ol 
“N de2 


dv 


01 


de 


+ V 


01 


The value of inductor current ij^-j is given by 


dv 


^X1 de 


01 


= Vj^lsin(e-eo)l 

'^01 

R. 


( 2 . 8 ) 


(2.9) 


The range of applicability of e is from to 0-| , corresponding to the beginning 
and the end of the conduction interval, so that 


ixi(0o^ = "'xi^®!^ ^ ^ (2.10) 

The solutions to (2.8) and (2.9), subject to the constraints given by 
(2.10), yield the following expressions for v^^ and i^^ : 

Voi(e) = Arsine + B^cose + exp(Be)[C^sin(a9) + D^cos(ae)] (2.11) 

ivi(0) = -^ {A.sine + B.cose + exp(3e)[C.sin(ae) + D.cos(ae)]}(2.12) 

A I Kj^ 11 11 

for 6 q £ e £ e-j . The values for the various coefficients depend on k, and 
6 q, and expressions for them are given in the glossary of symbols. 

When the diodes are not conducting, the equivalent circuit is that shown 
in Figure 2.3(B), where the second subscript 2 denotes this operating condition. 
The output voltage current this part of the cycle are described 

by the following equations: 

ix2(s) " 0 (2. I'D 

where the range of e is from e-j to (Sq+it), corresponding to the interval during 
which the inductor current is zero. Solving (2.13) yields the expression: 
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Vq 2(6) = VQ2(0-|)exp[-(9-0-j)/wCRj_] (2.15) 

From the requirements for a steady state in the discontinuous mode, and from 
the requirement of continuity of the output voltage across capacitor C, the 
following conditions must be satisfied: 

(2.16) 

(2.17) 

(2.18) 


’^01 




In order to obtain the steady-state solutions for Vq and i^^ the two 
transcendental equations given by (2.11) and (2.15) must be solved simul- 
taneously subject to the conditions imposed by (2.16), (2.17), and (2.18). 

The implicit relationships in this type of simultaneous equation pair suggest 
that digital computation be employed to obtain numerical answers for the 
solutions. An algorithm, described in the section which follows, was de- 
vised to rapidly determine the steady-state voltages and currents. 


2.2.3 Algorithm to Determine Steady-State Values of Vq and i^^ 

A simple algorithm for quickly determining the periodic steady-state 
is illustrated in the flow chart in Figure 2.4. Two first-trial initial 
values of output voltage and are entered as starting values 

and the corresponding initial trial values of 0g, and are 

calculated using the relationship 
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= arc sin (Vq^^VVj^) 

The values of 0^^^ and 0^^^ are then incremented by a prescribed small 

amount from the starting values of and respectively, and 

the corresponding values Vq^^^{ 0) and Vq^^^( 0) and 

are computed using (2.11) and (2.12). If the values of ix^^^ o'" 

are computed to be equal to or less than zero, the values of 0^^^ and 

0^^^ are identified tentatively as 0-j^^^ and 0-j^^^, which mark the 

end of the inductor conduction interval. As 0^^^ and 0^^^ are incremented 

over the interval from to + it and 0-j^^^ to + tt, 

respectively, and Vq^^^ are computed using (2.15). The inductor 

currents, sod ix^^^ remain zero during this interval. The error 

is computed as the difference between the computed value of 

at 0^^^ = ©Q^^^ + If and the initial guess for at 0^^^ = ©q^^^- 

Similarly, the error is computed for the second trial. Using 

and and initial values and a new trial value for 

fsl (3l f3l 

Vq^ ^ at O' ^ = ©Q ^ is computed by linear interpolation or extrapolation 
as shown in flow chart. The above process is repeated until the solution 
for VQ(0p) converges, thus yielding the values of 0 q and 0^. Results 
obtained from a number of computer runs show that usually only three 
or four iterations are required to obtain solutions which converge to with- 
in a prescribed error bound of O.OOOlVj^ . By incorporating this algorithm 
in the computer program, considerable reduction in the computation time 
involved in the steady-state analysis of rectifier-filter circuits when 
they are operated in the discontinuous-conduction modes, is obtained. 
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With these data the steady-state time waveforms of Vq and i^ can be 
calculated. The peak-to-peak output ripple voltage and the peak inductor 
current are recorded, and computer calculations are made for the average 
output voltage and the rms inductor current. Computer calculations are also 
made for the input power factor of the transformer-filter circuit. 


In the next section, an expression for the critical condition is 


derived. 


2.2.4 Critical Condition 


As mentioned in an earlier section, the critical condition establish- 
es the boundary between the two modes of operation. Using (2.2), the in- 
ductor current i^^ in the continuous mode is given by the expression: 

i -L cos(29-tl;) 

X ~ 7T 2R^ 3lZ(j2w)l 

where jZ(j2w)| and are the magnitude and phase of the filter impedance seen 
from the source vj at the second harmonic of the sinusoidal ac supply vol- 
tage and are defined in terms of k, Rj_, and coj^ as follows: ^ 




1 - + 4ic^ 


( 2 . 20 ) 


If) = arc tan [2k/ (1 - 4o3|^^)] - arc tan [Zio^j^/k] 


( 2 . 21 ) 


From (2.19), the minimum instantaneous value of i^ occurs when the 
cosine term has a value of unity. For the critical -condition boundary be- 
tween the continuous- and discontinuous-conduction modes of operation, the 
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minimum value of ij^ just reaches zero or 


4V 


Im ( 1 


= 0 


TT \ 2 \~ 3|Z(j2oj) 

Rearranging (2.22) and substituting for jZ(j2a))| from (2.20) leads to the 


( 2 . 22 ) 


relationship: 




2 4 2 

1 - + 16tdj^j + 4 k 2 

1 + (4u3j.jVk^) ^ 


(2.23) 


Squaring (2,23) and simplifying yields the expression: 

, X 4 2, R 2 4, 

F{k,(Ojj) = 0 = k + k (^ - 2tOj^ + 4(ii)j^^ ) g- 


(2,24) 


For a given value of the value of k satisfying the critical condition 
F{ic,(Dj^) = 0, referred to as k^^, can be found. Since both k and have a 
physical interpretation only for pos-itive values, for a given value of 
K is single valued, A plot of k for various values of ojm is shown in 

i w r ii 

Figure 2.5. Jhe abscissa, = toVLC, depends only on filter parameters and 
the frequency of the sinusoidal ac supply voltage. The value of k, pre- 
viously defined as wL/R^^, additionally depends on the filter' load resistance 
Rj^. For heavy loads, i.e,, Rj^ small, the value of k is greater than 
K^^, causing the filter circuit to operate in the continuous mode and the 
circuit is analyzed using Fourier-series techniques. For light loads, the 
value of K is less than k^^, the rectifier-filter circuit operates in the 
discontinuous mode, and the filter behavior is analyzed using piecewise- 
1 inear differential equations. The value of k^^ is approximately equal to 
it asymptotic value of 1/3 for values of greater than 5.0. The critical 
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Figure 2.5. Plot of Versus and Identification Regions. 
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condition and the circuit-analysis methods for continuous and discontinuous 
modes which have been described are used in a later section for generating 
design curves. 

2.2.5 Relationship Between ic, Q, ^ , and uj 

The transfer function of the filter circuit operating in the contin- 
uous mode, shown in Figure 2.2, is given by (2.3) and may be written in the 
familiar form Vq(s)/Vj(.s) = l/[(s/o3g)^ + (1 /Q)(s/cjq) + 1] where s is the com- 
plex frequency variable, wg = 1/ ^ LC is the familiar undamped resonant fre- 
quency of oscillation in radians per second, and 0 = R^/cogL is the quality 
factor at the resonant frequency. The damping factor ^ is related to Q by 
C = 1/2Q. From the previously given definition of W|^, the damping factor 
can be expressed as 

5 = k/2(i3|^ (2.25) 

The quality factor Q is related to k and 0 )^ by the expression: 

Q = V-K (2.26) 

The relationships shown in (2.25) and (2.26) should be helpful in identifying 
and attaching extended meanings to the design parameters < and 

2.2.6 Power Factor 

The input power factor for a nonlinear ac load such as either of the 
transformer-filter circuits of Figure 2.1 is defined by Schwarz [6], and by 
Kornrumpf and Walden [17] as the ratio of the power dissipated in the network 
to the product of the rms values of the ac input voltage and the ac input 
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current. Because of symmetry considerations and because the rectifiers and 

filter elements have been assumed lossless, the product of the ac input rms 

values (*^p/M 5 )VjgIpg equals the product of the rms values at the input to 

the LC filter Also, the output power to the load calculated as 

2 

Vge equal to the input power. Thus, the expression 


Power Factor = 


he he 


(2.27) 


has been used in calculating this performance measure of the network. 

2.3 Design and Performance Curves and Approximate Relationships 
2.3.1 Design and Performance Curves 

Using the methods of analysis described in the previous sections, a 
set of curves useful for both designing and analyzing rectifier LC-filter 
circuits was generated by digital computation and is shown in Figure 2.6. 

In the upper curves, the ratio of peak-to-peak output ripple voltage 
Vpp to dc output voltage Vg is plotted against k for values of lying be- 
tween 2.0 and 10.0. The dashed line in this plot represents the critical 
condition obtained from a solution of (2.24) and separates the plot into 
regions of continuous conduction and discontinuous conduction. In the second 
graph, two families of plots are presented for the same range of values for 
(Dj^ and K. They show the ratio of the dc output voltage Vg to the peak value 
of input source ''in, and the power factor P.F. of the ac input to the recti- 
fier-filter combination. These curves, Vg/Vj^ and P.F: versus k, actually 
consist of multiple plots of curves for the same values of as used in the 
top graph but appear to approach two single curves because of their weak de- 
pendency upon As seen from the Vg/Vj^ curve, for rectifier-filter cir- 
cuits operating in the discontinuous-conduction mode, a larger dc output volt- 


0.01 


0 . 001 * 

'‘°1 
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i 
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age is obtained than that in the continuous-conduction mode for a given value 
of input supply voltage Other useful design curves, the ratio of the 

rms value of current in the inductor to the average dc load current Iq, 
and the ratio of the peak value of inductor current to Iq are plotted 

as the third and fourth graphs of the same figure. Again, although appearing 
as single isolated curves these are actually the result of superimposed mul- 
tiple plots over the range of previously stated. 

Additional use may be made of the' graphs in Figure 2.6 to calculate 
the values of the rms currents in the transformer windings. In terms of the 


rms inductor current useful relationships are: 


Center-tapped secondary: 
Single-winding secondary: 
Either configuration: 


^Se “ ^XeV^ 

^Se "" ^Xe 

Ipe = (Ns/Np) Ixe 


(2.28) 

(2.29) 

(2.30) 


A study of the curves at the top of Figure 2.6 shows that various 

combinations of k and C 0 |^ may be chosen to meet a specified (^Qp/Vg) ripple 
requirement. Consequently, an additional specification must be provided to 
complete the information required for the design procedure. While the addi- 
tional constraint which usually is imposed is the desired output voltage, the 
peak or the rms value of the current in an individual diode, the maximum al- 
lowable current in the transformer windings, or the minimum desired power 
factor, may be used. Having chosen the additional constraint, the design 
process can continue, leading to values for the remaining unconstrained recti- 
fier-filter parameters which can be determined from the other plots in Figure 
2.6. To illustrate the design procedure, several examples are worked out in 
section 2.4. 
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2.3.2 Approximate Analytical Design Relationships 

The design curves for various performance parameters, shown in Figure 
2.6, can be approximated by analytical design relationships given in Table 
2.1. The values for the performance parameters computed from these approxi- 
mate design relationships are reasonably close to the computed values or 
those read from the design curves. The approximate design relationships can 
be automated on a digital computer or electronic pocket calculator to yield 
quickly designs for LC-filters. These relationships, therefore, provide an 
alternate method to the reading of data from the design curves. In Table 2.1, 
the analytical expressions are given for design relationships » 

and P.F.(ic). For the continuous mode, these 
approximate design relationships- are obtained from the exact relationships 
given in Appendix A under the assumption that is much larger than 1 and k. 

For the discontinuous mode, a good approximation to the (VQp/Vg) curve 
shapes in Figure 2.6 is obtained if they are considered to-be straight lines 
parallel to each other and starting at the continuous mode point for k= 1/3. 
One such approximation for fo'" the discontinuous mode is shown in 

Table 2.1, and values computed from this approximate relationship lie within 
ten percent of the values read from the curves. As mentioned earlier, the 
other performance parameters, ^Xm^^O’ weakly de- 

pendent on and hence are approximated by polynomials of second degree 
passing through the points on the curves for these latter parameters and are 
functions of k only. The values computed for these latter parameters from 
the approximations given in Table 2.1 lie within five percent of values read 
from the curves. In addition to the design relationships for the performance 
parameters, expressions for inverse relationships t^|>j(VQp/VQ,K) , ’ 




Table 2.1. Approximate design and performance relationships 
for LC-filter operating both in the continuous - 
and discontinuous-conduction modes.' 
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^Xe‘^^0^’ k(P.F.) are also given. In the continu- 

ous mode, Vg/Vjj^ is independent of k and hence there is no inverse relation- 
ship ><(Vq/Vjj^) for the continuous mode in Table 2.1. 

These approximate design relationships should be used carefully de- 
pending on the desired accuracy of the design. 

2.4 Design Examples 

To illustrate the usefulness of the curves presented in the preceding 
section, examples of two different filter designs now are discussed. These 
examples illustrate how two rather different application requirements can be 
satisfied. The International System (SI) of Units is used in the examples. 

Exampl e 1 . 

The values of the inductor and capacitor elements in an LC-filter cir- 
cuit with nominal output voltage Vq = 5V, nominal load current 1^ = 10 A, and 
ripple ratio found. The supply voltage to the recti- 

fiers is provided by the secondary windings of a center- tapped transformer, 
the primary of which is excited by an ac supply of 117vrms at 60 Hz. The • 
inductor-capacitor filter circuit is to operate in the continuous-conduction 
mode. 

Given that the ripple ratio VQp/Vg = 0.05 and that the filter is in the 
continuous-conduction mode, the design procedure begins with the V^p/Vg 
curves in Figure 2.6. Choosing the curve for a3j^j = 2.5 as an example, the rip- 
ple constraint under nominal load is met with a value' of ic = 6.0. Using the 
equivalent load resistance corresponding to nominal output voltage and 
load current, Rj^ = 0.5 ohm, the values of L and C may.be obtained using the 
relationships ic =wL/R|^, and = wyTc, to yield L = 7.96 mH and C = 5526yF. 
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Having established values for k, L and C, the performance of the 
filter network and the determination of important parameters useful in com- 
pleting the design of the remainder of the complete power supply may be ob- 
tained. For example, if the filter circuit always is to operate in the con- 
tinuous-conduction mode, the minimum vlaue of k is approximately 1/3 corres- 
ponding to a load resistance 18 times the nominal load resistance. Since 
Vq/Vj^ is very nearly constant over this range, the minimum load current be- 
fore the output voltage begins to rise is (Iq ^ominal'^^^^ 0.56A. If the 
maximum load current is, for example, 1.5 times the nominal Iq value of 10 A, 
the corresponding maximum value k reaches is 6.0 (15/10) = 9.0. 

Using o)j^ = 2.5 and the range of k from 1/3 to 9.0, the performance of 
the filter circuit is obtained from the various curves in Figure 2.6 as fol- 
1 ows : 

Ratio - Load Condition 



Minimum 

Nominal 

Maximum 

K = wL/R|^ 

1/3 

6.0 

9.0 

Average output voltage, Vq 

5.02 V 

5.02 V 

5.02 V 

Peak-to-peak ripple voltage to output 


• 

- 

voltage, 

0.055 

0.05 

0.043 

Output voltage to peak input voltage, Vq/Vj^^ 

0.64 

0.64 

0.64 

Average load current, Iq 

0.56A 

10. OA 

15.0A 

Peak reactor current to load current, 

2.0 

1.03 

1.0 

Rms reactor current to load current, I^e^^O 
Rms transformer secondary current to load 

1.2 

- 1.0 

1.0 

current, l 5 g/lQ 

0.85 

0.71 

0.71 

Power factor of filter circuit 

0.73 

0.90 

0.90 
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Of particular interest is the performance curve for the ratio of output volt- 
age to peak input voltage Vq/Vj^ which approaches the asymptotic value 
2/1T-0.64 for values of k greater than 

From the knowledge of the ratio l^^g/Ig and now 

possible to determine the characteristics required of the center-tapped trans- 
former and the diodes in the power supply circuit. Using the nominal load 
voltage and current values for this example, it is noted that the diode for- 
ward voltage drop {approximately 0.7 for silicon diodes) is a signifi- 

cant fraction of the output voltage Vg. Recognizing that in the continuous- 
conduction mode the average output voltage must be the average input voltage 
to the rectifier-filter less the constant diode drop, the required peak volt- 


age on the transformer secondaries must be ~ [Vg + For this 

examp 1 .95 V. Using the specified primary supply voltage, the turns 

ratio Nj/Np is computed to be 1/18.5. 

Using the nominal value of Ig as the maximum value for sustained oper- 
ation, the peak reactor current ly is approximately 10.3 A, while the rms 


reactor current 

Xe 


10.0 A. 


The rms secondary winding current in a cen- 


ter-tapped secondary is equal to the rms reactor current 1^^^ divided by 


and has a value of 7.07 A. The rms primary winding current, neglecting magne- 


tizing current and using a turns ratio of 1/18.5, is 0.54 A. The power factor 
P.F. is approximately 0.9, approaching closely the asymptotic value 2 x/F/tt. 


The diodes because of their shorter thermal time constant should be selected 


on the basis of the peak current to be encountered 10.3 A. 
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Example 2 . 

In this example, an LC filter circuit driven by a full -wave bridge 
rectifier connected directly to the ac supply is to operate in the discontin- 
uous-conduction mode. The specifications are: nominal output voltage Vq = 

115 V, nominal load current Iq = 5 A, with a ripple ratio Vgp/VQ = 0.02. The 
ac supply voltage is 117 V rms, 400 Hz. 

Under nominal output voltage and load conditions, the ratio Vq/Vj^ = 
0.695 and from the curve for this parameter in Figure 2.6, k = 0.2. To meet 
the ripple requirement Vqp/Vq = 0.02, the curve corresponding to (i)j^ = 3.5 is 
selected. Using the value of = 23 ohms for the nominal load, the values of 
the filter components are computed to be: L = 1.83 mH; C = 1060pF. 

If the filter circuit is to operate in the discontinuous-conduction 
mode for all loads within the range of the curves in Figure 2.6, k can vary 
from 0.01 to 1/3. Examining the characteristics of the combined transformer- 
rectifier-filter circuit over this range in k yields the following results: 

Ratio Load Condition 




Minimum 

Nominal 

Maximum 

K = (oL/Rj^ 


0.01 

0.2 

1/3 

Average output voltage, Vq 
P eak-to-peak ripple voltage to output 


152.2 V 

115.8 V 

105.8 V 

voltage, VQp/V^ 


0.0018 

0.02 

0.027 

Output voltage to peak input voltage, 


0.92 

0.70 

0.64 

Average load current, Iq 


0.33A 

5.0A 

7.6A 

Peak reactor current to load current. 

^Xn/^0 

4.6 

2.3 

2.0 

Rms reactor current to load current, 


1.9 

1.3 

1.2 


continued 
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Ratio Load Condition 

Minimum Nominal Maximum 

Rms ac supply line' current to load 

current, 1.9 ' 1.3 1.2 

Power factor of filter circuit, P.F. 0.68 0.74 0.73 

Using the nominal load voltage and current values as before, addition- 
al information on the characteristics of the remainder of the power supply 
circuit may be obtained. Since the diode forward voltage drop is small com- 
pared with the output voltage in this example it will be neglected. The peak 
reactor current and the peak current through the rectifier diodes is 

11.5 A (5.0 X 2.3). The rms reactor current is found to be 6.5 A (5. Ox 1.3) 
which is also the value of the rms ac supply line current for the full -wave 
bridge combination. The power factor P.F. is 0.74. 

Although the two examples which are provided to illustrate the use of 
the design and performance curves in Figure 2.6 are specified for either all 
continuous-conduction mode or all discontinuous-conduction mode operation 
over the load range, there is no reason for restricting the operation to one 
or the other modes. The designer has the choice of specifying the nominal 
load conditions and ripple requirement to fall in either mode and to enter 
the other mode under changing load conditions, provided the require- 

ment can be met with appropriate turns ratio. 

2.5 Conclusions 

r 

Methods for analyzing full-wave rectifier LC-filter circuits were pre- 
sented and these methods were used to develop design procedures for the rec- 
tifier-filter combination. The procedures cover the cases of operation in 
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both the continuous-conduction mode and the discontinuous-conduction mode. 

An algorithm was presented which permits rapid determination of the periodic 
steady state in the discontinuous mode using digital computation. Using the 
results for both modes of conduction, a set of design curves was generated, 
providing a straightforward procedure for designing both the inductor and 
the capacitor for a rectifier-filter circuit operating in the discontinuous 
mode and continuous mode. The use of the parameter K=toL/R|^ and the estab- 
lishment of a boundary condition on < for separating the two modes of opera- 
tion were shown to be useful in choosing operating conditions. 

From the design curves, currents in the different windings of the 
transformer and in each diode can be calculated. Using these curves, a de- 
signer can select the inductor, the capacitor, the diodes, and the trans- 
former to meet various applications requirements. Some of the requirements 
are: varying ac input supply magnitude, maximum allowable ripple, maximum 

rms and peak diode and transformer-winding currents, minimum pov^er factor 
and load conditions. Two examples were presented to illustrate the design 
procedures using the design curves. Results of experimental and simulation 
runs using the inductor and capacitor design values agree closely with the 
results predicted from the design procedures. 



CHAPTER III 


ANALYSIS AND DESIGN OF LC-FILTERS 
WITH NON-ZERO SOURCE IMPEDANCE 

3.1 Introduction 

In Chapter II, an analysis and design procedure for the LC-filter 
circuit operating both in the continuous and discontinuous modes, with the 
sinusoidal voltage source impedance assumed to be zero, is presented. In 
this chapter, the basic concept used in Chapter II is extended to include 
LC-filter circuits with non-zero source impedance. Figure 3.1 shows the 
transformer-rectifier-filter networks of Figure 2.1 connected to a sinusoidal 
voltage source through a primary-circuit impedance consisting of a series 
resistance and a series inductance |_j. This series resis- 

tance-inductance combination represents the Thevenin equivalent impedance 
of the usually encountered source and line impedance as seen from the trans- 
former primary-winding terminals. The transformers, shown in Figure 3.1, are 
assumed to be ideal. As mentioned earlier, the inclusion of source impedance 
in the design algorithm yields a more accurate and complete design. For 
each value of the dimensionless parameter p, defined q.s Rj/R[_ where Rj is 
the value of source resistance referred to the secondary side of the trans- 


( 34 ) 



Figure 3.1. LC-filter circuit with non-zero source impedance (^) R 

Np 2 

and (|^) Lj supplied by (A) center-tapped transformer 
rectifier and (B) full -wave bridge rectifier. 
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former, a set of design curves, similar to one shown in Figure 2,6 for the 
case with zero-source impedance, is obtained. 

The format of the presentation here follov/s closely that of the pre- - 
ceding chapter. The derivations of the useful expressions are presented in 
Appendices C and D. 

3.2 Analysis 

As in Chapter II, in the continuous-conduction mode, linear circuit 
analysis techniques are employed, while in the case of the discontinuous mode, 
the method of analysis requires computer solutions of the piecewise-1 inear 
differential equations. 

3.2.1 Continuous-Conduction Mode 

Assuming ideal diodes, the equivalent circuit for both of the confi- 
gurations in Figure 3.1 is given in Figure 3.2 along with typical waveforms 
of inductor current and output voltage. The two dimensionless parameters 
toj^] and K defined in Chapter II are redefined as w Jl^^C and where 

L = Lt + L is the sum of the source inductance referred to the secondary 
side of the transformer and the filter inductance. As described in Chapter 
II, the equivalent input voltage Vj may be approximated by the following 
expression when the circuit is operated in the continuous-conduction mode: 

v| = ( y - i C0S26) (-3.1) 

i IT CO 

Using the expression in (3.1) and the transfer function of the filter, 
which includes Rj and Lj, the useful relations for various voltage and cur- 
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(B) 


Figure 3.2. (A) Equivalent diagram for the rectifier-filter circuits 

with non-zero source impedance operating in the ; 

continuous-conduction mode. ' i 

I 

I 

(B) Associated voltage and current waveforms from com- J 

puter solutions of circuit equations. \ 

{ 

i 

I 


38 


rent components can be derived as in Chapter II with the results: 
DC output voltage: 


''o- 


Z''lm 

tt(p+1 ) 


(3.2) 


Peak-to-peak ripple voltage: 


8V 


Im 


■Op 


^ ^ (p+1 -4ca^) + { 2k+( Zpaj^/ic) ) ^ 


(3.3) 


Rms inductor current: 


4V 


Im 


Xe ttR, 




4(p+l)^ 18{(ptl-44)^ + (2 k+(2p4/k))^} 


(3.4) 


Peak inductor current: 


4V 


I 


Im 


Xm irR, 


1 


1 

21^ ■ 3 


1+(44/k^) 


(p+1-44)^ + (2k+(2p4/!c) )^ 


(3.5) 


3.2.2 Discontinuous-Conduction Mode 

The equivalent circuit when diodes are conducting is the same as for 
the continuous mode and is repeated in Figure 3.3(A). The output voltage 
Vqi is obtained from the solution of the following equation: 


2 ^^^01 


dv^ 


o VJ V ^ T 

4 + (p+Dvqi = Vjj^|sin(e-6Q)| (3.6) 

do 


The value of inductor current i^.| is given by 


^’xi de 


(3.7) 





Figure 3.3. Equivalent diagrams for the rectifier-fil ter circuits with 
non-zero source impedance operating in the discontinuous 
-conduction mode during interval (A) when diodes are con- 
ducting and (B) when diodes are not c;onducting. (C) 
Associated voltage and current waveforms from computer 
solutions of circuit equations. 
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2 

The solutions to (3.6) and (3.7) for the values of (k+(po)j^/k:)) 

2 

<4(jjj^j ( p+l ) , yield the following expressions for and ij^-|: 

Voi(e) = Ay^sine+By^cose + exp(g^e)[Cy^sin(a^e)+D^i„cos(a^6)] (3.8) 


ixi(0) = 1^ A^.^sine+B^.^cos0+exp(g^9){C.^sin(aj^0)+D.^cos(a^e)}^ (3.9) 

p p 

For the values of [ic+(pa}^/K:)] >4(,j^(p+l ) , the expressions for Vg-j and 
ij^.| are given by the following: 


Vgi(0) = Ay^sin 0 +B^^cos 0 +C^^exp(a-| 0 )+D^^exp(a 26 ) 


(3.10) 


''Xl^®^ R^_ 


A^ ^s ine+B^. ^cos9+C^I ^exp(a^ 


e)+D^l^ex 



(3.11) 


The values for the various coefficients depend on k, p, and 0g, 
and expressions for them are derived in Appendices C and D. 

When the diodes are not conducting, the equivalent circuit is that 
shown in Figure 3.3(B). The output voltage Vg 2 and current ij ^2 this part 
of the cycle are described by the equations (2.13) and (2.14) given in 
Chapter II. The expression for Vg 2 is obtained, as given 1n Chapter II, 

Vg2(0) = Vg2(0-|)exp[-(0-9^)/a)CRj_] (3.12) 

The steady-state solutions for Vg and i^ are obtained by solving 
the transcendental equations given by (3.8), (3.9), (3.10), (3,11), (3,12), 
and (2.14), as described in Chapter II. 
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3.2.3 Critical Condition 

The critical condition, k = corresponds to the boundary between 
the two modes of operation. Using (3.1), the inductor current i^^ in the 
continuous mode is given by the expression: 


4V 




Im 


1 cos(2e-i;; 

2R|_(p+l) ■ 3|Z(j2w) 


where |Z(o2o])| and ip are given by the following 


iZ(j2,,)| = 


{p+l-4o3^)^ + [2K+(2pco^/ic)]^ 
1 + (4(jj^/k:^) 


i/; = arc tan [2ic+(2pw^/K)]/(p+l-4a)^) - arc tan (2 co^/k) 


W 


(3.13) 


(3.14) 

(3.15) 


For the critical condition, the minimum value of ij^ just reaches 


zero or 


4V 


Im 


1 1 

2Rj_(p+iy ■ 3lz(j2<^) 


= 0 


(3.16) 


Rearranging (3.16) and substituting for lZ(j2(u)| from (3.14) leads to the 
relationship: 


/ (p+l-4w^)^ + [2 k+(2poj^/k)]^ \ 2 


(3.17) 


Squaring (3.17) and simplifying yields the expression 


F(ic,a)j^) = 0 = ^(p+1)^ - 2ai^+4tij^ - (4+8p-5p^) (3.18) 
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For a given value of o>j^, the value of k = which satisfies the critical 
condition F(k,o)j^j) = 0 given by (3.18), can be found. A plot of versus 
(Dj^ for various values of p is shown in Figure 3.4. For a given p, the value 
of approaches asymptotically ('^^+2p-l .25p^)/3 for values of greater 

than 5.0. 

3.3. Design and Performance Curves and Approximate Relationships 
3.3.1 Design and Performance Curves 

Using the methods of analysis described in the previous sections, sets 
of curves, similar to the one shown in Figure 2.6, for various assigned 
values of p = Rj/R^ weregenerated and are shown in Figure 3.5 through Figure 

3.12. Following the method in Chapter II, these curves may be used to cal- 
culate various performance factors for LC-fi Iters with non-zero source im- 
pedance. As seen from the design curves shown in Figure 3.5 through Figure 

3.12, for a particular value of k and the value of the ratio of peak-to- 

peak ripple voltage to average output voltage Vgp/Vg slowly increases with an in 

crease in the value for the parameter p = Rj/R^. The values for the other 
performance parameters P.F., and decrease with in- 
creasing p. The dependence of the parameters Vg/Vj^, and 

on p is greater in the discontinuous mode than in the continuous mode. The 
dependence of the other two parameters Vgp/Vg, and P.F. on p is greater in 
the continuous mode than in the discontinuous mode. In the continuous mode, 
the parameter Vg/Vj^^ is most sensitive to variations >n p and Ij^g/Ig is 
least sensitive. In the discontinuous mode, the parameter Vg/Vj^^ is again 
most sensitive to variations in p and P.F. is least sensitive. Thus, for 

particular values of the filter components L, C, and R^^, the ripple ratio 
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=o»/LeqC 


Figure 3.4. Plot of k versus u., for various 
^ cr N 

values of p = Rj/R|^. 
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cuLeq 


Figure 3.10 Design and performance curves for LC- filters with 
p = 0,3 operating both in the continuous- and dis- 
continuous-conduction modes. 
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V^p/Vo calculated from the design curves for the zero source impedance case 
shown in Figure 2.6 is smaller than that obtained from the design curves 
presented in this chapter. The values for Vq/Vj^ and P.F. read from the 
curves in Figure 2.6 are greater than thoseobtained from the design curves 
which take source impedance into account. Thus, the design curves pre- 
sented in Chapter II give optimistic values for the ripple ratio 
output voltage to peak input voltage ratio Vq/Vj^, and power factor P.F. 

For the filters with values of p other than those given in Figure 3.5 
through Figure 3.12, the design should be accomplished using one of the de- 
sign curves shown in Figure 3,.5 through Figure 3.12 corresponding to the 
value of p greater than and closest to the specified p. To illustrate the 
design procedure, an example is worked out in Section 3.4. 

3.3.2 Approximate Analytical Design Relationships 

For filters operating in the continuous-conduction mode, the design 

curves for various performance parameters, shown in Figure 3.5 through Figure 

3.12, can be approximated by analytical design relationships given in Table 

3.1. These, approximate design relationships are obtained from the exact 

2 

relationships given in Appendix C under the assumption that is much 
larger than 1 and k. In Table 3.1, analytical expressions are given for the 
design relationships Vqp/Vq(k,wj^,p) , Vq/Vj^^Cp), Ij^g/lQ(K,p) , ’ 

and P.F.(k,p). In addition to the design relationships for various perfor- 
mance parameters, expressions for the inverse relationships conj(Vgp/VQ,K,p) , 

Xm/^O’P) given. In the con- 

tinuous mode, Vg/Vj^ is independent of k and hence there is no inverse 
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relationship given in Table 3.1. The expression for another in- 

verse relationship k(P.F.,p) is also omitted because of the complexity of 
the expression. 

For the discontinuous case, the performance curves for each value of 
p can be approximated by the design relationships given in Table 3.2. The 
dependency of the various expressions on p is by means of coefficients a, b, 
a-| , a^, a^, b-j , h^, b^, Cp c^, d^ d2, and d^. The values for these 
coefficients for various values of p are given in Table 3.3. 

A good approximation to the (Vqp/Vq) curve shapes in Figure 3.5 
through Figure 3.12 is obtained if they are considered to be straight lines 
parallel to each other and starting at the continuous mode point for 
K = ( ^ l+2p-l .25p2 )/3 . As seen from Figure 3.5 through Figure 3.12, the 
other performance parameters, Vq/Vj^^, P.F., and are weakly 

dependent on and hence are approximated by polynomials of second degree 
passing through the points on the curves for these latter parameters as 
functions of k only. The values computed from the approximate relationship 
for (Vqp/Vq) lie within ten percent of the values read from the curves while 
the values computed from the approximate relationship for the other perfor- 
mance parameters, ^o'^^Im’ within five percent 

of the values read from the curves. 

3.4 Design Example 

To illustrate the usefulness of the performance curves presented in 
the preceding section,’ an example of a filter design is discussed. 



Table 3,2 Approximate design and performance relationships 
for LC-f liters with non-zero source impedance 
operating in the discontinuous-conduction mode. 


Vq 


T — U,p) 
^0 


P.F(ic.p) 






(dQ 


a-| - a^log <+ 32(109 


- b2logK+ b2(log 


- C2log'<-*- CjClog <)~ 


d-| - d^log K+ d2(log k)' 


10 

(logCc 



(®2 ‘ 

^a2^-Ka3C(VQ/Vjj^)-ai3 \ 

10 

\ 

2d3 } 

( 



\ 

TO 


2b2 1 

( 


.yj C2^->-4 c3[{ -c^ ] \ 

\ 

10 


2 C 3 / 


(h- 

.Jd2^+4d3C{P.p.)-d^3 

10 

\ 

2d3 / 


(V„/V, ,p) 






,p 


ParametsrsdQ&Q.a^ ,a2,a2,b2 ,b,,b^,c^ ,C2,C2,d^ ,d2,and d. 
are constants which depend on p = See Table 3.3. 
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3,4.1 Example 

The values of the inductor and capacitor elements in an LC~filter 
circuit with nominal output voltage Vg = 30 V, nominal load current 
Ig = 10 A, and a peak-to-peak ripple voltage to average output voltage 
of 4% are to be found. The supply voltage to the rectifiers is provided 
by the secondary windings of a center-tapped transformer, the primary of 
which is excited by an ac supply of 117 V rms at 400 Hz. The secondary 

C M \ ^ 

_2. I p "'s given 

as approximately 4.5 ohm and source inductance/ ] , is approximately 

\\) ‘•I 

80 yH. 


Using the specified primary and secondary supply voltages, the turns 
ratio Mg/Np is computed to be 1/3.9. The value of the source resistance 
referred to the load side Rj , is approximately 0.3 ohm. The value of source 
inductance referred to the load side Lj is 5.26 yH. Under nominal output 
voltage and load conditions, the ratio Vg/Vj^^ = 0.707 and from the set of 
performance curves for p = Rj/Rj_ = 0.1 in Figure 3.8, k = 0.09 from the 
Vg/Vim curve and the filter circuit operates in the discontinuous mode. 

To meet the nominal ripple requirement of Vgp/Vg = 0.04, the curve cor- 
responding to lOj^ = 2.0 is selected which provides a somewhat smaller value 
of ripple. Using the definitions of k = o^Lg^/R^ and = w and the 
value of R^ = 3 ohms for the nominal load, the values for the and C 
are computed to be: = 0.1074 mH, C = 5896 yF. The values of the filter 

components are computed to be: L = = 0.1022 mH; C = 5896 yF. 

If the filter circuit is to operate in the discontinuous-conduction 
mode for all loads within the range of the curves in Figure 3.8, k can vary 



from 0,01 to 0.363. Examining the characteristics of the combined trans- 
former-rectifier-filter circuit over this range in k yields the following 
results: 
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Load Condition 


Ratio (Using Non-Zero Source Impedance Curves) 

Minimum 

Nominal 

Maximum 


0.01 

0.09 

0.363 

Average output voltage, Vq 

31.8 V 

30.1 V 

24.6 V 

Peak-to-peak ripple voltage to, output voltage, 




’'op/''o 

0.0046 

0.032 

0.096 

Output voltage to peak input voltage, 

0.75 

0.71 

0.58 

Average load current, Ig 

1.2 A 

10.0 A 

33.1 A 

Peak reactor current to load current, 

3.3 

2.6 

2.0 

Rms reactor current to load current, Ixe'^^O 

1.6 

1.42 

1 .26 

Rms transformer secondary current to 




load current, Igg/Ig 

1.13 

1 .01 

0.89 

Power factor of' filter circuit, P.F. 

0.74 

0.77 

0.72 


Using the nominal load voltage and current values as before, addi- 
tional information on the characteristics of the remainder of the power 
supply circuit may be obtained. Using the values ~ 2.6 and 

(Ixe/^0^ = 1.42 read from the performance curves for k = 0.09, the peak 
reactor current and the peak current through the rectifier diodes is 
26 A and the rms reactor current is found to be 14.2 A. The rms secondary 
winding current in a center-tapped secondary is equal to the rms 

reactor current divided by JY and has a value of 10.0 A. The rms 

primary winding current, neglecting magnetizing current and using a turns 
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ratio of 1/3.9' is 3.64 A. The power factor P.F. is approximately 0.77. 

Since the diode forward voltage drop is small compared with the output 
voltage in this example it will be neglected. When the output voltage is 
small, however, the diode forv/ard voltage drop can be a significant fraction 
of the output voltage Vq. As mentioned earlier, the diode forward voltage 
drop was neglected in the analysis since this factor would have added an 
additional parameter to the family of performance curves presented earlier. 
There seems to be no easy way to take diode forward voltage drop into account 
for the case when the filter is operating in the discontinuous-conduction - 
mode without adding an additional normalized parameter to the design curves. 

The above example is now repeated for comparison purposes, using 
the design curves which assume zero-source resistance which were presented 
in Chapter II. 

Starting again with the ratio = 0.707, from the curve for this 

parameter in Figure 2.6, k is approximately 0.2. To meet the minimum ripple 
requirements ~ 0.04, the curve corresponding to = 3.0 is selected. 

Using the value of R|^ = 3 ohms for the nominal load, the values of the 
filter components are computed to be: L = 0.238 mH; C = 5986 yF. 

Examining the characteristic of the combined transformer-rectifier- 
filter circuit over the range of < in the discontinuous mode yields the 
following results 

Load Condition 


Ratio (Using Zero-Source Impedance Curves) 

Minimum 

Nominal 

Maximum 

K = (j)L/Rj^ 

0.01 

0.2 

1/3 

Average output voltage, 

39.1 V 

30.1 V 

27.2 V 

Peak-to-peak ripple voltage to output 




voltage, Vq /Vq 

0.0024 

0.028 

0.038 
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Load Condition 


Ratio (Using Zero-Source Impedance Curves) 

Minimum 

Nominal 

Maximum 

Output voltage to peak input voltage. 




yh,. 

0.92 

0.71 

0.64 

Average load current, Iq 

0.65 A 

10.0 A 

15.1 A 

Peak reactor current to load current, /I„ 

Xm 0 

4.6 

2.3 

2.0 

Rms reactor current to load current, 

1 .9 

1.34 

1.23 

Rms transformer secondary current to load 




current, 

1 .34 

0.95 

0.87 

Power factor of filter circuit 

0.68 

0..74 

0.73 


The peak reactor current and peak reactor current through the 
diodes is 23 A. The rms reactor current is found to be 13.4 A. The rms 
secondary winding current is 9.47 A. The rms primary winding current is 
similarly given as 3.44 A. The value of inductor L calculated using the 
design curves in Chapter II is greater than that needed in the design 
which takes source impedance into account, although the required value of 
capacitance C essentially remains the same. Thus the peak and rms values 
of reactor current, and the ripple ratio smaller for the filter 

designed using zero-source impedance curves than that designed using the 
curves presented in this chapter. For further comparison using the values 
for the filter components L and C as computed using the zero-source impedance 
curves, i.e. , k = 0.2 and = 3.0, the value for Vq/Vj^ from the curves 
for p = 0.1 in Figure 3.8 is found to be 0.64. Thus the actual dc output 
voltage Vq obtained using the filter designed with the curves presented in 
Chapter II instead of the curves for p = 0.1 is 27.1 V instead of the 
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desired 30 V. 

An example using the approximate analytical design relationships 
is presented next. 

3.4.2 Example Using Approximate Design Relationships 

The same filter design example is repeated here using the approximate 
analytical design relationships given in Tables 3.1, 3.2, and 3.3. 

Starting again with the ratio = 0.707, from the relationship 

for in Table 3.2, and using the values for a-j , agj and a^ 

corresponding to p = 0.1 from Table 3.3, k is computed to be 0.087. Using 
the specified value of ~ 0*04 in the relationship for “f^(VQp/Vg,ic,p) 

in Table 3.2, the value for is found to be 1.7. Using the value of 
= 3 ohms for the nominal load, the values for the and C are computed 
to be: Lgq = 0.1039 mH; C = 4404 pF. The values of the filter components 
are computed to be: L = 0.0983 mH; C = 4404 pF. 

Calculating the characteristic of the combined transformer-rectifier- 
filter circuit over the range of k in the discontinuous mode yields the 
following results: 

Load Condition 


Ratio {Using Approximate Analytical 
Design Relationships) 

Minimum 

Nomi nal 

Maximum 

= "-l-eq/RL 

0.01 

0.087 

0.363 

Average output voltage, Vg 
Peak-to-peak ripple voltage to output 

32.1 V 

29.9 V 

24.5 V 

voltage, Vgp/Vg 

0.0069 

0.039 

0.125 


Output voltage to peak input voltage, 


0.755 


0.706 


0.578 
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Load 

Condition 


Ratio (Using Approximate Analytical 
Design Relationships) 

Minimum 

Nominal 

Maximum 

Average load current, Iq 

1.22 A 

9.98 A 

34.1 A 

Peak reactor current to load current. 


• 



3.29 

2.64 

2.01 

Rms reactor current to load current, 




o 

1 — r 

tt) 

X 

1.62 

1 .43 

1 .25 

Rms transformer secondary current to 




load current, I 5 g/Ig 

1.14 

1 .01 

0.88 

Power factor of filter circuit, P.F. 

0.74 

0.76 

0.72 


The results obtained for the three examples for the nominal load 
conditions are summarized below. 

Example A: Using Zero Source Impedance Curves 

Example B: Using Non-Zero Source Impedance Curves 

Example C: Using Approximate Analytical Design Relationships 


Example 


Ratio 

A 

B 

C 

. = 

0.2 

0.09 

0.087 

Average output voltage, Vq 

Peak-to-peak ripple voltage to output voltage, 

30.1 V 

30.1 V 

29.9 V 

''op/''o 

0.028 

0.032 

0.039 

Output voltage to peak input voltage, 

0.71 

0.71 

0.706 

Average load current, Iq 

10.0 A 

10.0 A 

9.98 A 

Peak reactor current to load current, 

2.3 

2,6 

2.64 

Rms reactor current to load current, Iy„/Iri 

1.34 

1 .42 

1 .43 
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Exampl e 

Rati 0 A B C 

Rms transformer secondary current to load 

current, Ijg/lQ 0.95 1.01 1.01 

Power factor of filter circuit, P.F. 0.74 0.77 0.76 

It is seen from the above summary that the results obtained from 
the approximate analytical design relationships given in Table 3.2 are in 
good agreement with the results obtained using design curves presented in 
this chapter. 

3.5 Conclusions 

Methods for analyzing full-wave rectifier LC-filter circuits operating 
in the continuous and discontinuous modes with non-zero source impedance were 
presented in this chapter. The inclusion of source impedance adds an addi- 
tional normalized parameter p = Rj/R|_ to the set of design curves which were 
developed in Chapter II. Using the results of the analysis for both modes 
of conduction, a set of design curves for various assigned values of p was 
generated. 

4 

From the design curves, various performance factors of the filters 
can be calculated. The source impedance is always present in an actual 
physical system and its inclusionas a parameter yields a more complete 
design. An example was presented to illustrate the design procedures using 
the design curves. To illustrate the effect and importance of source im- 
pedance on the design and performance of LC-fi Iters, the same example was 
repeated using the zero-source impedance design curves presented in Chapter 
II. Using the values of sc and obtained from the design procedure using 



the curves for zero-source impedance, and working backward' through the 
design curves which take source impedance into account, shows that with 
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the values of k and obtained from the design curves in Chapter II, the 
dc output voltage ratio Vg/Vj^ will be smaller than the required value, 
although the ripple ratio Vgp/Vg will be smaller than that for the filter 
designed using source impedance curves because of the larger inductance 
value used. In general, as seen from the design and performance curves for 
the LC-filter circuits with non-zero source impedance, for specified values 
for the filter components L, C, and i.e., k and the ripple ratio 
Vgp/Vg is greater than that obtained from the design curves which do not 
take source impedance into account. The actual value of the parameter 
Vg/Vi^ is also smaller than that obtained from the design curves presented 
in Chapter II. Thus, a filter design obtained by use of the zero source 
impedance curves yields a set of parameters leading to an optimistic design 
in terms of ripple ratio and output voltage to input voltage ratio. In 
other words, higher values for the filter components L and C, and supply 
voltage are needed than those calculated from the zero-source impedance 
design curves to meet the specified requirements of ripple ratio Vgp/Vg 
and desired output voltage Vg. 



CHAPTER IV 


DESIGN OF AIR-GAPPED MAGNETIC-CORE INDUCTORS FOR 
SUPERIMPOSED DIRECT AND ALTERNATING CURRENTS 


4.1 Introduction 

Various procedures have been reported in the literature for the de- 
sign of air-gapped inductors carrying direct current. In most cases, a 
particular core with known size and shape is first selected for evaluation. 

The procedure then involves calculating a trial value for the number of 
turns and the length of the air-gap for the selected core and then analyzing 
the combination to determine if all the requirements, electrical and mechani- 
cal, for the inductor are met. If not met, then either the number of turns 
or the length of the air-gap is changed from its previous value and the 
calculations are repeated. The results again are checked against the require- 
ments, thus giving rise in most cases to an iterative trial -and-error proce- 
dure applied to the selected core. 

All design methods require the use of a considerable amount of mag- 
netic core data which usually can be obtained from core manufacturers. One 
of the earliest design methods is that presented by Hanna [16] which removes 
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the trial -and-error process applied to the selected core by using specially 
prepared magnetic core data. This procedure, however, is limited to induc- 
tors with small ac flux excursions and also requires the use of special 
magnetic core curves which are not readily available for some modern core 
materials. Two recently developed design procedures by Ray and Sartori [18] 
and by Mitchell [19] make use of the computer to eliminate the tedium usually 
associated with the normal trial -and-error process. An important and dis- 
tinguishing feature among the various methods is the way in which the core 
incremental permeability data are employed. The methods of [16] and [18] 
make use of a single incremental permeability curve which corresponds to the 
condition of very small ac flux excursion. The starting point in these pro- 
cedures is the specification of the desired inductance L and the direct cur- 
rent Ipj. in the winding. The method given in [19] and the one presented in 
this chapter make use of a family of incremental permeability curves and thus 


allow for large ac flux excursions. These procedures also start with a speci- 
fied value of inductance operating at a given value of direct current in the 
winding. However, additional design constraints consisting of the frequency 
and rms value of the fundamental component of ac voltage across the winding, 
f and respectively, are imposed on the design procedure. By making use 
of the known geometry of the selected air-gap type core, and the specified 
dc resistance and current-density specifications of the winding, the number 
of turns N are calculated in [18] and [19]. An optimum air gap corresponding 

4 * 

to maximum inductance of the core with winding of turns H is computed by 
the design procedure in [19]. The value of air gap found in [18] is 
not optimum in the sense defined above. The procedure by Hanna [16] is pre- 
sented in detail in the next chapter. 
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1 

J 

In this chapter two procedures for designing inductors are presented. 

I 

The first, working from a stored bank of magnetic-material data and core- 
geometry data, uses a computer to produce a design which is optimized in 
the sense of minimum core volume and minimum number of winding turns to meet 
a required minimum inductance over a specified range of circuit operating 
conditions. The second procedure is somewhat simpler; it does not yield a 
design which is optimized as defined above, but it has the advantage that 
the computations may readily be carried out on an electronic pocket calcula- 
tor. Both design procedures are applicable to most modern magnetic core 
materials including, for example, grain-oriented silicon steel, supermendur 
and permalloy. While type C cores are used in the examples presented in this 
paper, the basic procedures can be extended rather easily to other shapes 
such as type E cores, ring cores, and certain lamination structures. 

4.2 Design Requirements 

An inductor for filtering the ripple voltage in dc power supplies 
usually carries a direct current with an alternating current superimposed. 

In most cases, the requirements are that the value of the inductance over 
the anticipated range of circuit operating conditions be no smaller than 
some specified value which will provide for adequate filtering of the fun- 
damental component of a waveform of frequency f. For such a filter reactor, 
the ac flux excursion usually is large and is determined by the magnitude 
and frequency of the fundamental component of voltage-across the reactor. 

In the procedures developed in this paper, the following specifications are 
used for the design of filter reactors: 
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(1) L, the desired minimum value of inductance in henries. 

(2) the maximum value of direct current in the winding in amperes. 

(3) the minimum value of the effective or rms value of ac voltage 
at the fundamental frequency f across the winding in volts. 

(4) f, the frequency of the fundamental component of voltage appearing 
across the reactor in hertz. 

It can be shown that a filter reactor designed to provide an induc- 
tance L at a direct current and a fundamental ac voltage increases 
in inductance if either the direct current in the winding decreases or the 
ac voltage across the winding increases. Therefore, a design procedure for 
minimum specified inductance at a specified maximum direct current and mini- 
mum ac voltage will yield a reactor design that satisfies the specification 
requirements. 

4.3 Analysis 

In any design problem, a necessary first step toward finding a design 
procedure is an analytical study of the system to be designed. In the develop- 
ment that follows, the subscript m stands for quantities related to the 
magnetic-core material and the subscript g stands for quantities related to 
the air-gap. 

The magnetic circuit of the type C core shewn in Figure 4.1 is 
analyzed under the following assumptions. First, the length of the air-gap ■ 
Jig is quite small compared to the mean magnetic path length ji^ of the selected 
core; therefore, the total mean magnetic path length made up of Ji^ and ji^ is 
taken to be equal to i . Second, since also is small compared to the 
core build-up and the width of the magnetic tape, fringing flux is very small 
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compared to the total flux in the core; therefore, the effective cross- 
sectional area of the gap portion of the magnetic path is taken to be 
equal to the effective cross-sectional area of the magnetic material for 
the selected core. 

Having established the nature of the magnetic circuit, attention now 
is directed to the curves of Figure 4.2. Plotted in the figure are a nega- 
tive air-gap line relationship between and H, shown as Curve 1, and a 
normal magnetization curve relating and for a typical magnetic core 
material which is shown as Curve 2. The inductance of the reactor winding 
of N turns at a particular operating or quiescent point given by 


N^A. 


L = 


m 


(4.1) 


^0 

where Uq is the permeability of free space and is the relative incremental 
permeability of the magnetic material at the operating point [20]. An ex- 
pression for the minimum value of fundamental -component rms voltage 
across the reactor, is given by (4.2) in terms of the frequency of the fun- 
damental component f, the peak value of the fundamental component of the ac 

A 

flux density the number of turns N, and the cross-sectional area of 

the magnetic core A^^. 


/ac >/2'^''*^^m^mAC 


(4.2) 


These three relationships--the curves of Figure 4.2, Equation (4.1), 
and Equation (4.2)--can be manipulated, as shown in Appendix E, into an 
expression of the form 




(4.3a) 



-Curve2^Brn vs 


-Curve i: Bg=yU.o(NlDc--^rn 

Slope = -/ioim/^g 


50 


Magnetizing force H in A/m 


Figure 4.2. Curve 1; Negative air-gap line. 


Curve 2: Normal magnetization curve 
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The resulting expression derived as (ElO) is 



where the constants are defined as follows: 



(4.3b) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 


Ke 4 L/2yQA^ 


(4.9) 


In (4.3), which is the first of two basic relationships, all but the two param- 
eters and are obtained in a rather direct manner from the perform- 
ance specifications, from the quiescent operating point on the normal mag- 

4 

netization curve, and from the geometry of the selected core. 

The second basic relationship required for the design is obtained 
from the experimentally determined incremental permeability curves provided 
by the core manufacturer and is of the form ^ 2 (Bj^^Q,y^,H^g) = 0. The 
FzC*) = 0 relationship is graphically portrayed in Figure 4.3 where the re- 


lative incremental permeability y^ is plotted versus the peak ac flux density 
excursion for various values of the quiescent or bias-point values of 
H^q from the normal magnetization curve. Figure 4.3 is a plot of data taken 



Relative incremental permeability }±^ 


iOOOOf 


■ k-i^ 

- 0.47T 


^ = 79,577 





iOOO 


4 
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(K 


Peak ac flux 


Figure 4.3 Relative incremental permeal 
soidal flux density excursi( 
normal magnetization curve ■ 
steel . 
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from [21] and [22]. 

Using the analytical base just discussed, the framework for the design 
procedures is now developed. Assuming that a quiescent point 
a chosen core material and core geometry is assigned, solution of the two 
relationships Fi(*) = 0 and F2(*) = 0 is obtained simultaneously for 
and by iteration on a combination of curve reading and computation. Using 
these solution values, the number of turns N and the air-gap length 5,^ are 
determined, respectively, using (4.10), which is a rearrangement of (4.2), 
and using (4.11), which is derived in Appendix E as (E7). 

N = (4.10) 

t = 0.5K, (4.11) 

y 

where 

K7A-K1 - 4[K2 - (Kg/y^)] (4.12) 

The values of N and obtained for an arbitrary assignment of cor- 
respond to a design matching all of the imposed constraints except the. mini- 
mum turns requirements and a check for windability. To proceed toward an 
optimized design for minimum turns for a given set of operating specifica- 
tions, magnetic material and core geometry, an additional iterative operation 
is performed in which the quiescent operating point is varied until the solu- 
tion for N converges to a minimum value-. In a later section, two procedures 
for designing an inductor-one for an assigned quiescent point 
and the other for minimum inductor turns--are illustrated using a common set 
of specification requirements for the two designs. 
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4.4 Design Procedures 

The objective of both of the design procedures to be discussed is to 
determine the turns N, the air-gap length and the wire size which meets 
the designer's current-density specification for a selected magnetic core 
material and geometry. The algorithms used in these procedures are ex- 
plained in this section by reference to Flow Charts 1, 2, and 3 presented in 
Figures 4.4, 4.5 and 4.7. Individual parts of Flow Charts 1 and 3 will be 
identified by letters enclosed in brackets, e.g., [A], to which reference 
will be made in the section illustrating the design procedures by means of 
an example. 

In Flow Chart 1, Figure 4.4, starting with the assumption that a 

suitable core has been chosen and that a suitable quiescent point 

has been obtained in some manner, the algorithm for computing the number of 

turns N and the air-gap length i is demonstrated. An integral part of this 

y 

procedure is the simultaneous solution of the two relationships F^(.) = 0 
and F 2(0 = 0 to obtain values of and which are used in computing 

N and Jig. The algorithm used to compute N and developed in this chart 
is used as a sub-algorithm in the two remaining charts. 

The second chart. Figure 4.5, illustrates a procedure adaptable to 
computer-aided design techniques, in which a stored library of magnetic 
and geometric data on commercially available cores is used as the primary 
data source. This procedure begins by using data for the smallest core in 
the library and an initial quiescent point to calculate N and 2-^, 

using the algorithm in the first chart. The procedure continues by itera- 
tively selecting additional quiescent points and re-computing N and 2-^, 
until a minimum value for the number of turns N is found, along with a 
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physically-realizable air-gap length. If the core under examination is in- 

i 

capable of functioning to provide the specified requirements, it is discarded 
and the next larger core is examined. 

The third flow chart. Figure 4.7, is used to explain a design pro- 
cedure which does not produce a minimum-turns design, but incorporates a core 
screening function which enables a designer to choose suitable cores and to 
design an inductor for a specified quiescent point, using the computational 
capability of an electronic pocket calculator. The screening operation in- 
volves a relatively simple computation which indicates cores for which a 
physically-realizable air-gap is possible. 

4.4.1 Flow Chart 1 . 

Returning to the first flow chart in Figure 4.4, a closer study of 
the algorithm for determining N and is made. The operations in this chart 
begin with the assignment of the specification requirements of L, 1^^, 
and f and the choosing of a magnetic core material and geometry and a 
quiescent operating point These steps carry the procedure to an 

entry point designated by the number 1 enclosed by a circle. This entry 
point. Node 1, is used by the design procedures illustrated by the two re- 
maining flow charts, where the algorithm which follows Node 1 returns to 
these procedures values of N and £g. Since the values of and are 

assigned before entering Node 1, parameters Kj^ through Kg are calculated in 
Step [A] and treated as constants until another core and/or quiescent point 
is selected. The algorithm for finding N and depends on obtaining a solu- 

A 

tion set for and by an iterative evaluation of Fi(*) to produce a 

A, , 

computed error e by using trial values of and the corresponding values 
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of which are read from the incremental permeability curve for the se- 
lected value of H^q. 

To provide guidance on the magnitude of e representing a good engineer- 
ing solution, a range factor Ae is calculated as the magnitude of the dif- 
ference between extreme values, and e . in Step [B]. The defining 

niaX rrn n " 

relationships for and are shown in the flow chart. Using the range 
factor Ae and the precision factor of 0.001, the computation of the error 
factor e proceeds in a closed loop using Steps [C], [D], and [E] until the 

A 

solutions for and converge to within the permitted error bound. 

Once the solution set is found, the values of M and are computed in Step 
[F] and the design for the specified requirements of chosen core material . 

and geometry and the selected quiescent operating point is completed. i 

j 

4.4.2 Flow Chart 2 . I 

In the procedure illustrated in the flow chart in Figure 4,5, an : 

t 

optimization algorithm searches for the combination of quiescent point and 
air-gap length which yields a minimum value for the number of turns N on 
the reactor. Because no general mathematical model is available which will 
fit the wide variety of magnetic materials in use, purely analytical opti- 
mization techniques cannot be applied. However, by the use of an iterative 
■procedure which employs multiply-dimensioned polynomial fits to experi- 
mentally obtained data, an optimized quiescent operating point and air-gap 
length can be determined. .- . 

The algorithm in this flow chart operates on a base of data for com- 
mercially available air-gap cores such as the type C geometry. The data ’ 

base includes core dimensions and sample points taken from the normal mag- 
netization curve and from a family of relative incremental permeability 

I 

i 
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curves for the designated core material. For each sample point [H^( I ) I)] 
selected as a data-base point on the normal magnetization curve, there is a 
corresponding member of the family of relative incremental permeability 
curves which also is stored as a set of sample points , 

Starting with the smallest core in the data library and the first 
sample point [H^(l ) ,B^(1 )] on the normal magnetization curve, polynomial 
interpolation is used to approximate intermediate values of on the corre- 


sponding incremental permeability curve [H^{1 ) ,§^^^^(1 ,J) ,u^(l ,J)] when 
seeking the solution value in steps [C]-[D]-[E] of Flow Chart 1 the 
first time that this sub-program, marked by the circled number 1, is en- 
countered near the top of Flow Chart 2 in Figure 4.5'. The value of turns 
N and air-gap lengthily are thus determined for this first sample point. 

This procedure is repeated for each of the M sample points on the- normal 
magnetization curve. 

The particular sample-value quiescent point [H^(K) ,B^(K)] corres- 
ponding to the sma-llest calculated value for turns- N is used as the starting 
point for another interpolation routine. This interpolation procedure is 
used in a search for that particular intermediate value of quiescent point 
on the normal magnetization curve which produces the minimum number of 
turns N. The interval searched spans the range covered by this starting 
sample point and the two adjacent sample points. For each intermediate 
trial value of quiescent point polynomial interpo'l ation is used to 
approximate the corresponding value of on the normal magnetization curve 
when preparing to compute the design constants K-j to Kg in Step [A] of Flow 
Chart 1 when the sub-algorithm marked Node 1 is encountered near the center 
of the flow chart in Figure 4,5. The interpolating scheme for obtaining 


81 


the values of and and the value of y^ corresponding to the 

trial which are used in Steps [B]-[C]-[D]-[E] in Flow Chart 1, is 

described with reference to Figure 4.5. The interpolation routine makes use 

of Lagrange's interpolating polynomial of sixth degree [23]. The data points 

used in establishing the required polynomials correspond in all cases to the 

six points closest to the most recently established trial point. Points 

marked by crosses on the relative incremental permeability curves shown in 

Figure 4.6 represent the sample points stored as a part of the data base. 

The six relative incremental permeability curves corresponding to the six 

sample points on the normal magnetization curve which are closest to the 

trial quiescent point magnetizing force are used for interpolation. The 

values of y and y . for intermediate trial values of H ^ approxi- 
'^A,max A, mm mQ 

mated by interpolation using, respectively, the six data points, illustrated 


I in Figure 4.6 by circles surrounding the crosses, which correspond to the 


sample point values of y. for B and § on the six closest re- 

^ ^ ^A mAC,max mAC,mm 

lative incremental permeability curves. To obtain a value of y^ correspond- 

A. 

ing to the trial value of indicated in Figure 4.6 as trial 

needed in Step [D] of Flow Chart 1, two interpolations are used. The first 
interpolation provides values of y^ corresponding to B^^^ trial 


these six closest incremental permeability curves. The location of the 
points corresponding to these six values of y^ are shown as circles in 
Figure 4.6. These y^ values are then used as six data points in a second 

interpolation to find an approximate value of y^ corresponding to the trial 

/N 

quiescent point and trial* procedure of two interpolations 

is repeated until a solution value of is obtained in Steps [C]-[D]-[E] 

in Flow Chart 1. The values N and corresponding to this trial quiescent 


°mAC,min ^mAC, fricl 


°mAC,max 


A 

Peak ac.flux dsnsiry 


Ficure U.o niustration of asa of data-base ncints to calculate and u, 

iijFin .^,n2x 

for a trial value of lyin^ bst\veen H^(.<-1) and Also, illus- 

tration of the tv.'O iterations used to obtain a value of for an 

arbitrary trial H - and 3^,^ 

ml j , i.n c 1 
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permaability cata points stored as a 
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by internolation for the six values of 
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m 
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point are computed and N is compared with the lowest previously calculated 
number of turns. This trial quiescent value H^q is then incremented by a 
prescribed amount and the procedure described above is repeated. 

After the minimum value for N and its associated air-gap are found for ! 
a selected core, the maximum rms winding current is computed, a wire size 
is chosen, and the winding factor is computed. If the core is windable, 
the pertinent Information on core size number, minimum turns air-gap 

for the optimum design ^g^opt’ ^^isscent operating point (H^Q,opt’®mQ,opt^ ’ ‘ 

wire size and winding factor are printed. If the core is not windable, 
or if it is desired to obtain a design for a larger core, the next larger •; 
core size is selected and the procedure performs another set of computations \ 

i 

to obtain a possible design. j 

I 

4.4.3 Flow Chart 3 . i 

1 

The third flow chart, presented in Figure 4,7, shows-a-procedure • 

which may be employed to design an air-gapped magnetic-core inductor, using 
the computational capability of an electronic pocket calculator. This pro- . 

I 

cedure does not produce an optimized design with minimum turns for a given 
core and set of specifications. However, by means of a relatively simple 
screening process, the procedure does permit a designer to determine which 
cores among those available to him are capable of yielding designs with phys- 
ically realizable air-gaps. Cores which pass this screening process then 
may be considered further by using the algorithm for determining N and z 

9 i 

flow charted in Figure 4.4. 

The design procedure begins with the specification of the design \ 

requirements in Step [G].and the selection of a candidate magnetic core 

' ’ i 


1 

1 
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material, calculation of a parameter tiie selection of a core 

geometry and choice of quiescent operating point ^s indicated in 

Step [H]. The screening rule, which is developed in Appendix F, depends 
on decisions made on inequality expressions involving a constant Kg which 
is a function of certain of the design specifications and of the core geo- 

I 

metry. In Step [J], if Kg is greater than the design equations 

will yield a workable core and N and are computed as indicated by the 
branching to Node 1, which represents the entry point to the N and 
algorithm of Flow Chart 1, If this Kg inequality relation proves to be 
false, then Kg is tested in Step [K] to determine if it is greater than 
Pg{H2Qy^Q)/2, where y^q is the maximum value of on the relative incre- 
mental permeability curve corresponding to the selected quiescent point. 

A true condition of this inequality again leads to the computation of N 
and Ag. The false condition leads to a final check on the value of Kg 
in Step [L] where it is compared to the parameter Po(*^^qU^q)fj,-jn^^* 
latter parameter is unique to the particular magnetic material chosen and 
is the minimum value of Uq(H^q^^q)/2 over the range of and y^g for 
j which data are available. If the value of Kg is greater than the value of 
I parameter nQ(^mqU^q)fg-in^2, the core will permit a design with a physically 
realizable air-gap if a new quiescent or bias point is chosen in Step [N] 
which satisfies the inequality in Step [K]. If the inequality in the third 
decision inequality is not satisfied, then a larger core must be selected 
in Step [M] and the process repeated. 

When a combination of core geometry and quiescent operating point is 

found which leads to values of N and 5 , , the rms winding current is deter- 

y 

mined in Step [0] and a wire size is chosen according to the current density 
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constraint imposed by the designer. The winding factor then is computed in 
Step [P] and if the winding factor is acceptable, Step [Q], the design is 
completed. If the core is not windable, a larger core is selected and the 
screening process is re-entered. 

4.5 Example to Illustrate Design Procedures 

To illustrate the design procedures previously described, a sample 
design for a filter inductor is presented' in outline form. The procedure 
which is illustrated in some detail is that shown in Flow Charts 1 and 3, 
Figures 4.4 and 4.7, and is suited to the use of an electronic pocket cal- 
culator on which the computations shown were made. 

Flow Chart 3 . 

[G] Specifications: L = 0.5 H; 1^^^ = 0.15 A; = 1.0 V; f = 60 Hz; 

*^w,max ” 

[H] (i) Choose magnetic material : 

Grain-oriented 3 percent silicon steel 

Normal magnetization curve (Figure 4,2, Curve 2) 

Relative incremental permeability curves (Figure 4.3) 

Calculate (H^nu n) • ^ Fi^om extreme left end of upper curve 
' mQ^AQ min 

in Figure 4.3, (0.1 x 79.577)2(2800) = 1 .77 x 105(A/m) 

(ii) Select Core: 

= 0.120015 m; A^ = 1.613 x 10'‘»m2; A^^ = 4.838 x 10~‘^m2 
(iii) Select a bias point: 

A bias point is arbitrarily assigned by the designer. Choice 
of a bias point in the knee region of the normal magnetization 
curve usually has been found to lead to workable designs. Choose 
= 31.83 A/m (0.4 oersted); = 1.2 T (12.0 kilogausses) 


1 

1 


87 

i 

WmQ - 


[I] 

Calculate K„: K = 2.9072 x 10^ 

o o 


[0] 

Kg > go to Node 1 , the entry 

point for the algorithm for 


calculating N and 2-g. Since Kg > steps [K] and [L] are not 


used. 


Flow 

Chart 1 


t — 1 
< 

1 1 

Calculate K^, K 2 .' Kg, K^, Kg, Kg: 
Ki = -5.2863 X 10'® 

Kg = 1.6003 X 10'^^ 


Kg = 7.3047 X 10"^ 

K 4 = 5.409 X 102 


Kg = 0.24003 

Kg = 1.2334 X 10^ 

[B] 

Read from relative incremental permeability curve. 


For = 31.83 A/m (0.4 oersteds) 

= 6200 

a ,max 

at B np = 1.4 T 

mAC ,max 


V, • = 700 

A, mm 

at B^.p =-0.02 T 

mAC jim n 


Calculate: e . , Ae 

max mm 

= 2.1908 X 10^ 
max 

= -6-039 X 10^ 
mm 

i 

Ae = 8.2305 x 10^ 


1 — I 

0 

1 1 

A 

Choose trial value of 



A 


Trial 

No. 

^mAC T 

Next Step 

1 

0.04 

Go to [D], trial no. 1 

2 

0.03 

Go to [D], trial no. 2 

3 

0.0279 

Go to [D], trial no. 3 

4 

0.02807 

Go to [D], trial no. 4 


I 
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[D] Read value of from relative incremental permeability curve cor- 

A 

responding to H^q = 31,83 A/m and value of from [C]; calculate 
Ky and e: 


Trial 

No. 

% ^7 

e 

Next Step 



1 

880 2.4242 x 10“*^ 

2.9737 X 105 

Go to [E], trial 

no. 

1 

2 

800 2.5098 X lO'^^ 

7.8649 X 10^* 

Go to [E], trial 

no. 

2 

3 

785 2.5275 x lO"'* 

-5.9982 X 103 

Go to [E], trial 

no. 

3 

4 

790 2.5215 X lO"'* 

-7.3500 X 102 

Go to [E], trial 

no. 

4 

UJ 
1 ■■■> 

Test value of jej against precision range- 

factor 



j 

O.OOlAe = 8.2305 x 102; 


» 



Trial 

No. 



Next Step 



t 

! 1 
1 

Isl = 2.9737 X 105 

> 8.2305 X 102 

Go to [C], trial 

no. 

2 

1 

2 

lei. = 7.8649 x 10^^ 

> 8.2305 X 102 

Go to [C], trial 

no. 

3 

3 

le| = 5.9982 x 103 

> 8.2305 X 102 

Go to [C], trial 

no. 

4 

i 

1 4 

t 

[el = 7.3500 X 1Q2 

< 8.2305 X 102 

Go to [F] 



\ 

Calculate N and 





1 

N = 828 turns 





1 

1 

[ 

z = 1.261 X 10“*+ m 
9 





1 Return 

to Flow Chart 3, step [0] 





i 

' Flow Chart 3 





[0] 

Calculate and choose wire size: 





= 0.150211 A 






Using a current density of 

5.0671 X 10"7 m 

2/A (1000 circular mils/A) 

9 


wire size AWG 28 is satisfactory. 



i 

i 
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[P] Calculate windabil ity: 

For AWG 28 with heavy Formvar insulation, the total area of the wire 
and insulation = 1.051 x 10 ^ (207 circular mils); F^ = 0.18. 

Selected core is windable and design is completed. 

The design example illustrated is not an optimized design, in the sense 
of minimization of turns. Using the same magnetic core material data and 

I 

core dimensions, the following results were obtained from a computer program 

which implements the optimization algorithm described earlier in Flow Chart 

2 of Figure 4.5. The results of the nonoptimized design results are recapit- ' 

1 

ulated for comparison. 

Optimized Design Nonoptimized Design | 

I 

N 738 828 ' 

I 

1 .156 X lO"'^ m 1 .261 x lO"'* m 

30.24 A/m (0.38 oersted) 31.83 A/m (0.4 oersted) • 

I 

1 

4.6 Conclusions i 

i 

The design procedures presented in this chapter are particularly useful j 
in the case of filter inductors for electronic power supplies where the wind- 
ing usually must carry a large direct current in addition to presenting a 
high reactance to the fundamental component of the alternating current. By 
incorporating into the procedures experimentally determined relative incre- 
mental permeability curves, the large ac flux excursions normally encountered 
in such applications can be taken into account. In the computer-aided de- 
sign procedure, an important feature is the optimization algorithm for deter- , 

I 

mining the quiescent operating point and air-gap length yielding the minimum 

i 

number of turns. An easily used screening process which permits a designer ! 

i 

i 
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to identify candidate cores for which a physically realizable air-gap is 
possible and a relatively simple procedure for solving for the incremental 
permeability and ac flux density excursion using an electronic pocket cal- 
culator are provided in the second procedure. 



i 
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CHAPTER V 


DESIGN OF AIR-GAPPED INDUCTORS USING SPECIALLY 
PREPARED MAGNETIC CORE DATA 


5.1 Introduction 

In Chapter IV it was mentioned that Hanna's method provides a direct 
and easily used procedure for the design of air-gapped inductors by using 
specially prepared magnetic core data [16]. Hanna's method, however, is 
limited to the design of inductors with small ac flux excursions. The meth- 
od uses specially prepared magnetic core curves which unfortunately are not 
available for some of the modern core materials. 

This chapter presents a direct calculation procedure for the design 
of air-gapped inductors with large ac flux excursions. The heart of the 
new design procedure depends on the development of specially prepared mag- 
netic core data calculated from core material data which are available 
from core manufacturers. Using the results of the analysis in Chapter IV, 
and following the pattern of Hanna's development, a set of design curves 
is generated. These curves make use of a family of relative incremental 
permeability curves and thus allow large ac flux excursions to be taken in- 
to account. For a selected core, by using these curves a designer can find 
the length of the air gap 5.^ and the number of turns N for the winding to 
meet the design requirements given in Chapter IV. 


( 91 ) 
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Using the results of the analysis of a magnetic circuit v/ith an air 
gap presented in Chapter IV, the design relationships are developed. Fol- 
lowing a discussion of the generation of design curves with the help of a 
digital computer, an example is given to illustrate the design procedure. 
The procedure for generating the design curves is presented in detail in 
flow chart form in Appendix H. Toward the end of the chapter a comparison 
is made between Hanna's curves and the design curves presented lathis 
chapter for corresponding cases of small values of ac flux excursions. 


5.2 Design Relationships Used in Generating Hanna's Curves 

The design relationships used in Hanna's method, given by (5) and 
(6) in Ref. [16], are repeated here as (5.1) and (5.2) 


,,2 8 ^.[— + X 10 "® 

LI PC mQ Upp 


V 


0.4ir(--^ + ^) 


(5.1) 


^ ^m 


Nipip Bp, , a 

DC _ niQ r 1 ^ _£-j 

a^ 0T4Tr‘-Ur,p a„*^ 

m DC m 


(5.2) 


where \) = a A is defined as the magnetic core volume. When these equations 
mm ^ 

are written in SI units, using "^^hey take on the form: 




V 


1 K 


(5.3) 


NIpp a^ ®mQ^ 


a_ 


m 


Ur 


(5.4) 
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Rearranging (5.3) and (5.4) leads to 

r2 


LI 


u 

V ^0 


1 + u 


/g % |. II ,2 
q m 0 


(5.5) 


A Z 


m 


MI nr 

DC _ g tnQ . n 
m 111*^0 ^ 


(5.6) 


The values for the relative incremental permeability y^ used in (5.5) by 

Hanna are for very small values of ac flux excursions and thus depend only 

on the dc operating point on the normal magnetization curve. A typical 

curve of y^ versus for 3-percent grain-oriented silicon steel for very 

small ac flux excursions, obtained from References [21] and [22], is shown in 

Figure 5.1. Using data from such a curve and from the associated normal 

magnetization curve, such as shown in Figure 4.2, Chapter IV, for an assigned 

value of air-gap ratio equations (5.5) and (5.6) are treated as para- 

2 

metric equations with as the parameter, to obtain a curve of LIqq/v 
versus MIpq/£jjj. For a given value of air-gap ratio, several values of H^q 
are assigned and the corresponding values of y^ and B^^q are used in (5.5) 
and (5.6) to determine the associated values of Ll^q/v and MIqq/) 1|^. A plot 
of these two quantities for various values of air-gap ratio is known as a 
set of Hanna's curves, and is illustrated for 3-percent grain-oriented steel. 
Permalloy 80, and Supermendur in Figures 5.2, 5.3, and 5.4, respectively. 

Using Hanna's curves, an inductor can be designed to meet the re- 
quirements of specified inductance L, and dc current Ipj. in the winding. For 

2 

the selected core, the value of LIpq/v is determined using the specified 
values of L and Ipq, and the known volume v of the chosen core. The corres- 
ponding value of air-gap ratio S, /i and NI^^/Jl are determined from the 

y 111 m 



Relative Incremental permeability fju^ 


Fi gure 



i 


5.1 . 


Relative incremental permeability versus dc operating point 
H„o on normal magnetization curve for small values of ac 
flux excursions for 3-percent grain-oriented silicon steel. 


t 




i 

' i 
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Figure 5,4. Example of Hanna's curves generated by computer program for 
Supermendur 



98 


set of Hanna's curves. The required value of air-gap I and turns N are thus 
determined from the known path length of the selected core, and the 
specified value of 

In the next section, Hanna's method is extended and design relation- 
ships which take large finite ac flux excursions into account are developed. 


5.3 Design Relationships when AC Flux Excursions 'are Large 


In this section, a set of design relationships for use in generating" 
design curves is developed. From the analysis of a magnetic circuit with 
an air gap presented in Chapter IV, the expression for N is given by (E-2) 
in Appendix E and is repeated here for convenience. 


. _ ^q^mQ Vo''’mQ 
■ ''O^DC 


(5.7) 


Substituting the value of N given by (5.7) into the expression for L given 
by (4.1) in Chapter IV, the resulting expression becomes 


>- = U, 




0 

m 


^O^DC 


)‘ 


(5.8) 
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Rearranging (5.8) leads to 
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(5.9) 


A 


'OC 


m 


Recalling the definition of the magnetic core volume A; = in (5.9) and 
rearranging leads to the first basic relationship given by (5.10). 

r 2 


LI 


DC 


= V, 


0 


i + H 

i lin 
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(5.10) 
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The second basic relationship is obtained from the expression for 
the alternating voltage across the inductor given by (4.2) in Chapter IV, 

= ^2irNfA„i^C (5.11) 

Rearranging (5.11) leads to 


^mAC 



VZTTNfA^ 


(5.12) 


The other basic relationships required for the design are obtained from the 
experimentally determined normal magnetization curve and incremental permea- 
bility curves as illustrated, for example, for 3-percent grain-oriented 
silicon steel in Figure 4.2 and Figure 4.3, respectively. These four rela- 


tionships, which include equations (5.10) and (5.12) and curves typified by 

Li 2 V.p 

Figures 4.2 and 4.3, contain seven variables DC., H » 

'* ‘‘ni ,/5iTNf/V, 

®mAC’ ^A' manipulated into a single relationship in 


terms of the four variables 


V ’ 


, B and which can be symbolical- 


ly represented by the functional notation 



V 



(5.13) 


More particularly, using symbolic functional notation, the value of used 
in (5.10) is obtained from the family of relative incremental permeability 
curves where the dependence of upon both the dc operating point on 

/N 

the normal magnetization curve and upon the peak ac flux excursion can 

yN 

be indicated by use of the symbol u^(Hj^q, Similarly, the value of 

B^g obtained from the normal magnetization curve as a function of H^^g only 
can be expressed functionally as expression for LI^^/v 

given by (5.10) can be written in a form useful as a design relationship 
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similar to the one given by (5.13) as 


V ^0 ^ Un mQ^ 




0 


(5.14) 


LI2 '•• 

which relates the four parameters ^C, B nn and H 

— 5^,^ mAC mQ 

A second relationship used in the design procedure to be discussed 
in the next section is obtained by rearranging (5.7) in the form 


NInp Jl„ B . 

DC = _a. jnQ. + H 

S'- 5-_ Pn 


m 


''m >^0 


(5.15) 


This relationship given by (5.15) is used in computing a design parameter 

NIpc 

(- ^" ■) ^ for the design curves. 

^m P A 

For very small values of peak ac flux excursions the design re- 

lationship (5.14) used in the present development reduces to (5.5) used in 
Hanna's method and the second relationship given by (5.15) is identical to 
the corresponding (5.6). The design curves which make use of design rela- 
tionships (5.14) and (5.15) and which take large finite ac flux excursions 
into account are presented next. 

5.4 Design Curves 


Families of curves for three different materials, 3-percent grain- 
oriented silicon steel. Permalloy 80, and Supermendur, which relate design 
parameters and for various values of air-gap ratio vmre 
generated with the help of a digital computer and are shown in Figures 5.5, 
5.6, and 5.7. These families of curves were obtained using magnetic mater- 
ial data in the form of the normal magnetization curve and the family of 














lidc 



Figure 5.7. Design curves for Supermendur which t 


finite ac flux excursions into account. 


o 

CO 
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relative incremental permeability curves in conjunction with the design re- 
lationships (5.14) and (5.15). The algorithms used in generating these 
design curves are explained in detail in Appendix H. For any assigned value 
of air-gap ratio S- the right-hand side of (5.14) can be considered to 
be a function of and only, and a curve of LI^^/v versus for 
some value of can be produced. To do this, for a specified a value 

/\. A 

®mAC ^mAC min assigned, and the values of and corresponding 

/N 

to H and B . are obtained from curves such as those of Figure 4.2 and 
mQ mAC,min 

Figure 4.. 3, respectively. These values are substituted in equation (5.14) 

2 

to determine the corresponding value of LIp^/v. Similarly, the values for 
7 

LIpc/v corresponding to for various values of are obtained in 

this manner. The optimum value of H^^g = H^g for a specified value of 

2 

a /% which yields the maximum value of LI^p/v is used in generating the de- 
g m uu 

sign curve for a given ^g/\- The value of = ^^^DC'^^m^opt 

ponding to H^^g^^p^ is calculated from (5.15). Using H^g^^ for H^g in 
(5.14) and varying B^^^g from the specified minimum B^g to the maximum 

A rt /N 

B nn > a curve of LIn«/v versus is obtained as illustrated in Figure 

mAC,max DC mAC 

5.5 for silicon steel. In a similar manner, design curves for other a-ir- 
gap ratios are obtained. 

Using these design curves, an inductor can be designed to meet the 
requirements of specified inductance L, dc winding current Igg, rms value 
of the sinusoidal inductor voltage V^g, and frequency f of the ac voltage. 

O 

Both the ordinate tlgg/v and the parameter (*^I[)Q/^tj,)opt ^ increase as the 

air-gap ratio is increased, as seen from the design curves. For a selected 

2 

air-gap ratio, the value for parameter Ll^g/v also increases with increase 

/-w 

in value of peak ac flux excursion For a selected core, the value of 



LIpj>/v is determined using the specified values of L and and the known 
volume V of the chosen core. The corresponding values of Z /Z^y ^^^DC^^m^oot’ 
and are determined from the curves in Figures 5.5, 5.6 or 5.7. The 
required values of air gap Z^ and turns N are thus determined from the 
specified values of Igj. and the known mean path length of the selected 
core. The actual value for the is calculated from (5.12). To meet the 
design requirements under all operating conditions, the value of B^^^ cal- 
culated from (5,12) should be greater than the value obtained from the de- 
sign curves; in other words, the designer should insure operation on the 

selected curve at a point corresponding to a value always higher than 

2 

the calculated LIp^/v. 

In the following section, the complete procedure for designing an 
air-gapped inductor using these design curves is developed and an example 
is presented. 


5.5 Design Procedure 


Using the design curves developed in the proceeding section, a 
step-by-step procedure for designing air-gapped inductors is now presented, 

[A] The specifications for the inductor are: desired minimum inductance 

L; maximum dc winding current Ip^; minimum rms value of ac inductor 
voltage frequency f of ac voltage; maximum winding factor 


w,max‘ 


[B] 1) Choose magnetic material 

2) Select core: v, window area A^^^. 

[C] Calculate LI^p/v: Using the specified value of L and Kp in [A] and 
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[D] 


CE] 


[F] 


known v of the selected core in [B], the value of LI^^/v is calcu- 
lated. 

Choose a design curve corresponding to the smallest air-gap ratio 

that intersects the line parallel to the abscissa and passing 

2 

through the ordinate point corresponding to LI^^/v calculated in [C]. 
Read from the design curve: Read the values of corresponding 

A 

and at the intersection. 

/N 

Calculate N, and From the known specifications and core 

dimensions in [A] and [B] and from [E], the air-gap length turns 

N, and B^^j, are calculated using the equations 

I = {£ /£ )l 
g ' g' m m 

(Nlnr/^ ) + ^ 

^ DC m^opt m 


I 


DC 


'AC 


"mAC 


V 27 rMfA, 


m 


[G] If calculated in step [F] is greater than that obtained in step 
[E], continue to step [H]. If not, choose a curve corresponding to a 
larger air-gap ratio and return to step [E]. 

[H] Determine wire size and check windability: Calculate the rms wind- 

ing current 1^^^ using step [0] in flow chart 3, Figure 4.7 in Chapter 
IV and select wire size based on this from a wire table. Compute the 


winding factor F 


, If F < F ^ , the design is complete. If not, 
w w = w,max ^ 


select a larger core and return to step [B], 


To illustrate the design procedure, an example is now presented 
using these design curves. To provide a comprehensive comparison of all 
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three methods, the same example is repeated using Hanna's curves, and using 
the computer program for optimum design presented in Chapter IV. 


5.5.1 Example Using New Design Curves 


[A] Specifications: desired minimum inductance L = 0.55 H; maximum dc 

winding current Ip^ = 0.2 A; minimum rms value of ac inductor volt- 


[B] 


[C] 

[D] 


'AC 
factor F, 


w,max 


age V^p = 5.0 V; frequency of ac voltage f = 120 Hz; maximum winding 
= 0.4. 

1) Choose magnetic material: 3-percent grain-oriented silicon 

steel 

2) Select core: % = 0.12 m; A = 1.613 x 10 ^ m^; v = 1 .936 x 10 ^ 

* .mm 

m^; A = 4.838 x 10 ^ m^. 
wn 

Calculate LIp^/v: LIqq/v = 1137 

The design curve corresponding to ^^ 1 ^^ ~ 0.0014 is chosen in Figure 

5.5. 


[E] Read from design curves for the chosen magnetic material: read the 

values for ^mAC Figure 5.5 corresponding to 

the value of LI^p/v =1137 calculated in step [C]. 

= 0.0014 
^m 


NI 


DC, 


( ) . = 1489 ampere-turns/meter 

op U 

m 


[F] Calculate N, and 

From the known specifications and core dimensions in [A] and [B] and 
from [E] 
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a = 1.68 X 10"^ m 

g 

N = 894 


^mAC 


[G] calculated in step [F] is greater than that obtained in step 
[E], continue to step [H]. 

[H] Calculate and choose wire size: From step [0] in flow chart 3 

in Chapter IV 


I„ = 0.2007 A 
Xe 

Using a current density of 5.0671 x 10 ^ m^/A (1000 circular mi Is/ A), 

wire size AWG 27 is satisfactory. For AWG 27, the total' area of the 

-7 9 

wire and insulation A^^ = 1.312 x 10~ m (259 circular mils); 

F^ = 0.24; selected core is windable and design is completed. 


5.5.2 Example Using Hanna's Curves 


The example in Section 5.5.1 using the same core size and magnetic 
material is now repeated using Hanna's Curves. Corresponding to the value 
of LIp(yv= 1137 calculated in step [C], the values for and are 

read from Hanna's curves in Figure 5.2. 

£ 

^ = 0.0016 


MIqC 

= 1550 ampere-turns/meter 

x»,., 

m 

The values for the air gap £g and N are calculated from the known 


core dimensions and specified Iqq. 
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N = 930 

ilg = 1.92 X 10”^ m 

Based on the dc current in the winding = 0.2 A and from step [F] 
and [6], the wire size AW6 27 chosen and winding factor is calculated as 
0.25. 

Thus it is seen that designs obtained using Hanna's curves require 
a larger air gap and larger number of turns N for the winding to meet the 
same design requirements. 

The above example was repeated using the computer program presented 
in Chapter IV. Using the same magnetic core material data and core dimen- 
sions, the following results were obtained from a computer program which 
implements the optimization algorithm described earlier in Flow Chart 2 of 
Figure 4.5. The results of the designs using Hanna's curves and the design 
curves presented in this chapter are recapitulated for comparison 

N g 


Hanna's Curves 

930 

1 .92 X 10"^ m 

New Design Curves 

894 

1 .68 X 10"^ m 

Optimized Design, Chapter- IV 

861 

1.60 X 10*^ m 


From the above table it is seen that the design obtained using 
Hanna's curves yields the largest number of turns N and air gap £g and the 
design obtained using the computer-implemented optimization algorithm yields 
the smallest number of turns N and air gap to meet the same design re- 
quirements. The design obtained using curves presented in this chapter re- 
quires a larger air gap and larger number of turns N than the optimized 
design but smaller than the design using Hanna's curves. 
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5.6 Comparison with Hanna's Method 

The design relationships used in generating Hanna's curves and the 
design curves presented in this chapter, are given by (5.5) and (5.6), and 
(5.14) and (5.15), respectively. As mentioned earlier, the value of rela- 


tive incremental permeability in (5.5) is for very small values of peak 

ac flux excursion and depends only on the dc operating point on 

the normal magnetization curve. The values of y^ used in (5.14) are taken 

from a family of relative incremental permeability curves shown in Figure 

4.3, and depend on both the snd H^g. Thus for a selected air-gap 

ratio, a curve of LIq^/v versus is obtained when the dc operating 

point Hj^g on_the normal magnetization curve is varied from some minimum 

value to some maximum value. As seen from Hanna's curves in Figure 5.2, for 

2 

a given air-gap ratio the curve of approaches a maximum value for 

some H^g. In the case of the design relationship given by (5.14), as ex- 

2 

plained in Appendix H, for a given air-gap ratio a curve of Ll^g/v versus 

/s 

B^Ac obtained for each dc operating point H^g on the normal magneti- 

zation curve as the peak ac flux excursion B^^g is varied from some minimum 

2 

to maximum value. For a given H^g, the minimum value of LIgg/v occurs at 

the point where B^^g is minimum. Hanna's curves are for very small values 

of ac flux excursions, and thus for a given' air-gap ratio, Hanna's curves 

2 

can be considered as loci of these minimum values of LIpg/v. As discussed 

2 

in Appendix H, only the curve corresponding to the maximum value of LIgg/v 

is used as one of the design curves presented in Fig^ 5.5 in this chapter. 

2 

For a given air-gap ratio, the minimum value of i-Igg/v that is obtained 
from the design curves in this chapter corresponds to the maximum value of 



in 


2 

LIdc/"^ obtained from Hanna's curves. Thus, for very small values of ac flux ! 

excursions, the design curves presented here approach the maximum values of 

Hanna's curves. By taking larger ac flux excursions into account, advantage 

is taken of higher values of relative incremental permeability in calculat- 

2 

ing the design parameter LI^^/v. Therefore, for the same value of desired 
inductance, a smaller number of turns is always required when using the de- I 

sign curves presented in this chapter than that determined using Hanna's ' 

curves. ' 

I 

5.7 Experimental Verification 

A number of type "C" cores made of grain-oriented 3-percent sili- 
con steel, have been designed and assembled according to the results from \ 

the design procedure. The measured values of inductance and ac flux- | 

density excursion agreed with values obtained from the design procedure ; 

I 

within five percent. Because it is impractical to adjust the air-gap ; 

I 

length to the exact design value specified by the design results, a part i 
of the difference between, the measured and predicted performance can be 
attributed to this effect. 

5.8 Conclusions 

1 

Specially prepared magnetic core data in the form of curves useful for • 

i 

designing air-gapped inductors were presented. These curves make use of i 

i 

r 

core data available from manufacturers of air-gapped cores. Making use of 
these curves, a straightforward method of determining the length of the 

r 

t 

- . \ 
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air gap and the number of turns for the winding is presented. As opposed to 
Hanna's curves, by incorporating into the design curves relative incremental 
permeability data, the large ac flux excursions encountered in filter induc- 
tors are taken into account. Advantage is taken of the larger values of 
relative incremental permeability for finite ac flux excursions. For a 
specified air-gap ratio, these design curves yield the maximum value for 
LIp^/v. The design obtained from Hanna's curves is more conservative as it 
yields a larger number of turns N for the winding than needed to operate 
under finite ac flux excursions. The design obtained from the design 
curves presented in this chapter is not optimized in the sense defined in 
Chapter IV but always yields a smaller number of turns than that obtained 
from Hanna's curves. 


1 

\ 



I 


CHAPTER VI 


DESIGN OF ENERGY-STORAGE REACTORS FOR DC-TO-DC CONVERTERS 
USING AIR-GAPPED MAGNETIC-CORE STRUCTURES 


6.1 Introduction 

This chapter develops the procedure for designing the energy-storage 
reactor for dc-to-dc converters using an air-gapped core. The development 
here follows closely the material presented in Ref. [15]. 

Methods for designing the energy-storage reactor for dc-to-dc con- 
verters using computer-aided and table-aided techniques have been developed 
which use energy balance relationships and a designer-selected maximum flux 
density in the magnetic core as the major design constraints [9, 12]. These 
procedures cover the reactor design for the twelve converter configurations 
which result from the combination of the three single-winding power stages 
(voltage stepup, current stepup, voltage stepup/current stepup), and' the 
two-winding voltage stepup/current stepup power stage, and the three con- 
trollers (constant frequency, constant on-time and constant off-time). Sim- 
plified circuit diagrams of the four power stages are illustrated in Figure 
6.1. Developed particularly with gapless powdered permalloy and ferrite 
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Figure 6.1. Four Energy-Storage DC-to-DC Converters. (A) Single-Winding Voltage Stepup. 

(B) Single-Winding Current Stepup. (C) Single-Winding Voltage Stepup/Current 
Stepup. (D) Two-Winding Voltage Stepup/Current Stepup. 
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toroidal cores in mind, past procedures make use of a collection of data 
from a finite population of core sizes and permeabilities. This paper 
builds upon the analytically derived design equations presented in these 
references and develops new procedures which permit the design of the 
energy-storage reactor using a magnetic core with an air gap. The use of 
such cores permits the selection to be made from an innumerable population 
of core size-permeability combinations and makes available a wider range 
of magnetic core sizes, particularly in the larger volume sizes. 

6.2 Analysis of the Magnetic Circuit 

The introduction of an air gap of length a into a magnetic circuit 

9 

with mean magnetic path length causes the permeability parameters usually 
associated only with the non-air-gapped circuit to change rather considerably. 
The composite-circuit permeability parameters, usually referred to as ef- 
fective permeabilities, may be obtained from approximations given below, 

I provided that the length of the air gap is small enough so that the total 
► 

path length of the air-gapped circuit, a + £ , may be approximated by the 

Ml y 

quantity and that thev cross-sectional areas of the magnetic material 
A and the air gap A are taken to be equal, thus implying negligible fring- 

I *** y 

ing in the gap region. Under these circumstances, the air gap modifies the 
effect of the relative normal permeability of the magnetic material, desig- 
nated by the symbol to produce an effective relative normal permeability 
for the composite magnetic circuit defined by 

^DC,eff " ^^DC^g^^m^^ 


(6.1) 


r 
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and the effective relative incremental permeability defined by 


"a.eff ° 

where is the relative incremental permeability of the magnetic material 
alone at the dc bias point of the core. If the quantities ^nd 

(yA^g/itj^) are both large with respect to unity, then these two effective, 
i relative permeabilities are equal and can be represented by the single sym- 
bol Ugffs where 

^eff " ^DC,eff " ^A,eff " V^g 


Since at any bias point the incremental permeability of the magnetic material 
is less than the permeability obtained from its normal magnetization curve 
at the same point, the condition necessary for (6.3) to hold is that 
I ^^A^g^^m^ large with respect to unity, or 


>> 


m g 


(6.4) 


This condition is one that normally is met without difficulty for 
most air-gapped magnetic circuits when the magnetic material operates below 
the knee of its saturation curve and the air-gap ratio (s-g/^,^) is moderately 
large. An important implication of this condition is that essentially all 
of the energy stored in the magnetic core is concentrated in the air gap 
and very little in the magnetic material. The air-gap.ped reactor appears 
as a linear element and thus permits use of the analytic relationships of 
Ref. [12], which are based on an energy-balance analysis of dc-to-dc con- 
verters, to be applied directly to the design of air-gapped reactors. 
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In Ref. [12], Eq. (19), a fundamental relationship between the re- 
quired minimum volume of a homogeneous powdered permalloy or ferrite magnet- 
ic core is established. This relationship, rewritten as 

'’min " ^O^'eff® 


is in terms of Pq, the permeability of free space, the effective relative 
permeability, and a factor s which is defined as 

«'24W^,n,ax/(Bn,ax-8R)^ . (6-6 

The quantity is the amount of energy in joules which must be trans- 

ferred by the reactor over one switching cycle under the operating condition 
corresponding to maximum power flow through the converter. This quantity 
i is a function of the power stage specifications and a time parameter v/hich 
depends on the type of controller. Parameter B is a designer-specified 

IIIQ A 

maximum value of flux density which is reached but is never exceeded. The 
choice of a value for depends on the type of magnetic material selected 
for the reactor. Bj^ is the residual flux density of the magnetic material. 

Using the expression given in (6.3) for effective relative permeabi- 
lity in an air-gapped core subject to the restriction set by (6.4), the 
minimum volume for meeting the energy transfer requirement for a converter 
is 


'’min " ^ 

Since the volume of an air-gapped core is equal to the actual cross-sectional 
area of the magnetic material (gross core area times stacking factor) times 
the mean magnetic path length, i.e., in order to meet the require- 
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merits of (6.7), the inequality 


Vg ■ V 

must be satisfied. This relationship is of fundamental importance in the pres- 
ent development. The left hand side, under the assumption = Ag, rep- 
resents the volume of the air gap portion of the composite magnetic struc- 
ture. Thus, (6.8) establishes a minimum value for the volume of the air gap 
in terms of the permeability of free space and the parameter 6 which is a 
function only of the power stage and controller specifications and the maxi- 
mum allowable flux density in the core. 

From the air-gap volume inequality of (6.8), it is seen that the de- 
signer has the choice either of picking a value for A^, using magnetic core 
catalog data, then computing the required minimum value of air gap £g, or 
picking a value for 2 ,^, then obtaining an appropriate core with A^ greater 

than the minimum value for core cross-sectional area A These two 

m,min 

choices are expressed in (6.9) and (6.10), where the use of minimum values 
allows an equality relationship. 


^g,min = ^0^/^m 
^m,min "" ^O^^^g 


(6.9) 

( 6 . 10 ) 


In summary, in this section the conditions required of an air-gapped 
core to function as the energy-storage reactor in a da-to-dc converter have 
been established which provide a simple relationship between the effective 
relative permeability of the composite structure and the mean magnetic path 
length and the air gap length. Using the results of previous procedures 
based on maximum energy transfer requirements, a relationship betv/een the 
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volume of the air gap portion of the core and converter specifications and 
maximum allowable flux density was developed. This relationship may be 
interpreted to yield minimum values for either the core cross-sectional 
area or the air gap length. In the following section, the complete pro- 
cedure for designing an energy-storage reactor for this application is de- 
veloped and an example is presented. 

6.3 Design Procedure 

Given the specifications of the power stage and controller, the 
essential elements of the design procedure involve the selection of the 
geometrical form for the magnetic circuit, the type of magnetic material, 
and the. determination of the core size, the air gap length, and the number 
of turns. In this section, the detailed steps are outlined and an example 
is presented to illustrate the design procedure. The power stages included 
in the design procedures are the three frequently encountered single- 
winding configurations for voltage stepup (VU), for current stepup (CU), 
and for voltage stepup/current stepup (UD), and the two-winding configura- 
tion for voltage stepup/current stepup (2UD). Fo,r each power stage, three 
controllers — constant-frequency (FQ), constant on-time (TN), and constant 
off-time (TF) — are considered. In Ref. [12], equations are identified, 
where appropriate, by equation numbers with suffix letters which identify 
the power stage type. The same form of identification is used when needed 
in the step-by-step procedure which follows. 

A. The specifications for the converters considered in this paper are: 

desired output voltage V^; expected minimum and maximum input voltages 

'^I min ^I max’ 't^'^^’^sistor saturation voltage V^; diode voltage drop V^; 

expected maximum value of output power P„ ; 

u * rn8 X 
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and a time parameter which depends on the type of controller employed — the 
switching period T = 1/f for (FQ), transistor on-time t^^^ for (TN), or 
transistor off-time t^^^ for (TF). In the two-winding configuration, the 
presence of the additional winding, as pointed out in [10] , provides an addi- 
tional degree of freedom in the design procedure. By utilizing the extra 
degree of freedom, any one of various design constraint options can be in- 
corporated into the procedure. These design constraints are: A, duty cycle 

centered at a specified value; B, specified minimum duty cycle; C, specified 
range of duty-cycle variation; D, specified maximum transistor collector- 
to-emitter voltage; E, specified maximum diode reverse voltage; H, specified 
maximum duty cycle; and J, specified turns ratio. Table II in Ref, [12] pre- 
sents the option constraint equations for the turns ratio Y " Ty “ ^p/N^, 
where the subscript y represents one of the constraints A, B, C, D, E, H 
or J. 

B. Choices that are to be made which are associated with the magnetic cir- 
cuit are: type of magnetic structure and material, and maximum allowable 

flux density Compute the values of ? and <s, which are functions of 

the converter specifications and using from Ref. [12], Eq. (1), (3), or 
(5), depending on the controller type (VU), (CU), (UD), or (2UD). 

C. Specify either the air gap length or the cross-sectional area A^ of 
a selected magnetic core. If the air gap is specified, continue to step D. 

If A^ is chosen, then additionally selecting the core with the smallest value 
of £ results in the smallest core volume and gives the Greatest assurance 
that the inequality in (6.4) is satisfied. Proceed to step E. 
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D. Compute the value of using the value of chosen, select a core 

with A greater than or equal to A . and a value of £ as discussed in 
m m,min m 

step C, and compute the value of Kg using the equations 

m,min 0 g 

Kq = 1 + [1 - (?/5)(A^ 

9 m,min m 

E. Compute the value of £^ using the value of A^ chosen, select an 

air gap £^ which is greater than or equal to £^ and compute the value 

of K„ using the equations 
y 

g,min ^0 m 

Kg - 1 + [1 - (./S)«lg_„.„/;lg)]'''2 

Continue to step F. 

F. Using the value of £ obtained in step C or E, and the value of £ for 

( 9 m 

the core selected in step C or D, compute the effective relative permeability, 
using the approximation given by (6.3) 

^eff " V^g 

As a caution, at this point the designer should be aware that too large a 

value for might indicate that the inequality in (6.4) is not satisfied. 

In such a case, choice of a larger value for £ or a smaller value for A^ 

9 , m 

in step C will result in a smaller Continue to step 6. 
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G. Compute K^^using the appropriate expression, according to type of power 
stage. 






(V2Po,max^(^x-^R) 


(VU) 

(CU) 

(UD) 

(2UD) 


f (V2Po,max^^Vl,min-V("max-^R)/(Vl,min"V''D-''Q^ 
^(V2Po,max^^'^I,min-'^Q)(Vx-BR)/^^(^,min-''Q^^ 

Using the value of obtained in step C or E, and the values of K^q and 

Ko just obtained, compute the turns N (for single-winding reactors) or N and 
y p 

Ng (for a two-winding reactor) and inductance L or and L^. 


N or N = K K £ /p„ 
P 9 10 g 0 



L = A 

^0 m g 


P 0 p m g 


4 = ^?W 


Compute the maximum rms winding current Iw„ using Ref. [12], 

Ac 5 iildX 

Eq. (2), (4), (6), depending on controller type, for (VU), (CU), (UD), or 
(2UD), Select the wire size and compute the windability of the reactor. 
If the winding factor is too large, return to step C and choose a smaller 
air gap or choose another core of same with larger value of £|^. 
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To illustrate the design procedure outlined above, an example design 
is now presented using the International System of Units (SI). The power 
stage is a single-winding voltage stepup configuration and is controlled by 
a constant-frequency controller. 

A. The specifications are: 


Umin = IS '' 

V„ = 0.5 V 

Q 

I .max = 2'* '' 

V. = 0.8 V 

D 

= 28 V 

f = 2 kHz 

n m,v ■ AOO W 
0 ,max 

T = 500 pS 


B. Choose a C-core structure with magnetic material of grain-oriented 3% 
silicon steel with tape thickness of 1.016 x 10”^ m (4 mils). Choose 
= 1-0 T. Using Eq. (1,VU), Ref. [12], assuming the residual flux 

TTlaX 

density Bj^ = 0, 


= 0.1543 


-4 

C. Choose the air gap length to be 6.096 x 10 m (24 mils). 

D. Compute A . , select core, and compute K„. 

m,min 2 

A^ “ 3.181 X 10”'^ m^ 

m,min 0 g 

-4 2 

Choose a core with gross core area of 3.632 x 10 m . Using a stacking 

-4 2 

factor of 0.9, A = 3.269 x 10 m . Choose the smalTest available value 
m 

-3 2 

of a = 0.1832 m and window area A,,„ = 1.116 x 10 m . 
m wn 

Kg = 1 . [1 - = 1.164 
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F. Compute 

y = J, /& = 300 

■^ef f m g 

G. Compute K^„,N, L, and I„ . For the voltage stepup configuration, 

iU A0 jITIqX 

So- (V2Po,max)(Vi,,in “ '^qH^max " ' ^q) 

= 2.164 X 10"^ 

N = SSoV^O^ 

>- = "oNV^g = 

Using Eq. (2,VU), Ref. [2], 


Xe,max 

’’o.n,axt''o"''D-''o> 


1 + 


12 


‘-f'o.n.ax'VVV' 

x-7 „.2 


1/2 


= 24.6 A 


Using a current carrying capacity of '5.067 x 10~ m /A (1000 circular mils/A) 
a wire size of AW6 6 is selected. For double-coated insulation, the winding 
factor is 0.167, which is satisfactory, and the design is complete. 


6.4 Conclusions 

The procedure for designing the energy-storage reactor for dc-to-dc 
converters which is reported in this chapter is applicable to a variety of 
types of magnetic structures in which an air gap is inserted. The analysis 
which leads- to the fundamental relationship for the required minimum volume 
of the air gap region depends upon the principal assumption that the product 
of the incremental permeability at the core operating point and the air gap 
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ratio be large with respect to unity. This condition usually can be met 
with most common structures and materials. Since the required minimum air 
gap volume is the product of the gap length and cross-sectional area, the 
designer has the choice of specifying either one and computing the minimum 
value of the other. The quantity determining the minimum required volume 
is a function of the power stage and controller specifications and the 
specified maximum flux density in the core. 

Previously reported research results provide the basic expressions 
which permit calculation of the minimum air gap volume for twelve combina- 
tions of power stages and controllers. Once the air gap volume is determined, 
either by initially specifying the air gap and choosing a core with given 
cross-sectional area or by specifying a core at the outset with a given 
cross-sectional area and then determining the air gap, the remaining important 
design parameters can be obtained. Using the values for the physical dimen- 

1 sions of the chosen core, including air gap, the effective relative permeabil- 
ity , turns, inductance, and maximum rms winding current can be determined, 
a wire size selected, and windability checked. If the chosen core is not 
windable, the design procedure leads back to a point where another choice 

t 

may be made which will lead to a windable design. 

Since the procedure is independent of the type of air-gapped structure 
and magnetic material and depends only on the physical dimensions of the 
chosen core and air gap, it should be useful in many applications where 
powdered permalloy or ferrite cores may not be applicable or available be- 
cause of costs or size requirements. 


CHAPTER VII 

CONCLUSIONS AND SUGGESTIONS 
FOR FUTURE RESEARCH 


7.1 Conclusions 

The research reported in this dissertation provides the power elec- 
tronics system designer with useful analytical guidance and practical usable 
procedures for designing LC-filters for use with full-wave rectifiers, and 
air-gapped inductors for filters and for a group of widely used dc-to-dc 
converters. The design procedures presented in this dissertation eliminate 
most of the conventional trial -and-error procedures which are frequently 
followed in designing these power subsystems. Methods for analysis of 
full-wave rectifier LC-filter circuits operating both in the continuous- 
conduction mode and discontinuous-conduction mode are also presented. The 

I design procedures presented here are well suited to computer-aided design 

I 

I 

techniques, but calculations required in most of the procedures are easily 
made on an electronic pocket calculator. 

Chapters II and III present the steady-state analyses of rectifier 
LC-filters operating both in the continuous- and discontinuous-conduction 
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modes. In Chapter II analyses are presented for the case when the impedance 
of the sinusoidal supply source is assumed to be zero. • In Chapter III 
a finite value for the source impedance is assumed. In the' continuous- 
conduction mode, linear circuit analysis techniques are employed, while 
in the case of the discontinuous mode, the method of analysis requires 
computer solution of the piecewise-1 inear differential equations. An 
algorithm is presented which permits rapid determination of the periodic 
steady-state in the discontinuous mode using digital computation. From 
these methods of analysis sets of curves and approximate relationships 
useful for designing filters are generated. Using these curves or relation- 
ships, a designer can select the inductor, the capacitor, the diodes, and 
the transformer to'meet a set of specified requirements. 

In Chapter IV, two design procedures are presented for the design of 
inductors using air-gapped magnetic cores. The first uses a computer to 
produce a design which is optimized in the sense of minimum core volume and 
minimum number of turns to meet a required minimum inductance over a speci- 
fied range of circuit operating conditions. The second procedure does not 

i 

! yield a design which is optimized in the above sense, but it has the advan- 
tage that it may readily be followed by hand calculations or with a calcula- 
tor. An easily used screening process which permits a designer to identify 
candidate cores for which a physically realizable air gap is possible, is 
also provided in the second procedure. 

Using the analysis of a magnetic circuit with air gap which was pre- 
sented in Chapter IV, specially prepared magnetic core data in form of design 
curves are generated and presented in Chapter V. The most important feature 
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of this approach is the straight-forward method for computing the number of 
turns and the length of the air gap for the selected core for meeting a set 
of design requirements. This approach does not yield as much information 
as the methods in Chapter IV do, but it does provide the designer with an 
easy method to quickly reach a workable design. 

In Chapter VI, a procedure for designing energy-storage reactors for 
dc-to-dc converters using air-gapped cores Is presented. The analysis, based 
on the energy-transfer requirement of the reactor, leads to a simple rela- 
tionship for the required minimum volume of the air gap. Determination of 
this minimum air gap volume then permits the selection of either an air gap 
or a cross-sectional area. Having picked one parameter, the minimum value 
of the other immediately leads to selection of the physical magnetic struc- 
ture , Although the magnetic cores used in the examples throughout the 
i dissertation are type C cores, the design equations and procedures presented 
are independent of core geometry and magnetic material used and can be ex- 
tended rather easily to other shapes such as type E cores, ring cores, and 
certain lamination structures. It is believed that the design procedures 
I presented in this dissertation- can significantly reduce the time and the 

I 

effort required to design LC-fi Iters for full -wave rectifier supplies, filter 
inductors, and energy-storage reactors in dc-to-dc converters. 

7.2 Suggestions for Future Research 

Three particular areas for future research hav& suggested themselves 
during the course of this research effort. The first area is concerned 
with the analysis of the LC-filter when it is a part of the regulator system 
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1 in which the source waveform is a square wave with a certain duty ratio. It 
is suggested that the analysis methods presented here be investigated for 
applicability to the complete filter-regulator system. 

Secondly, in the dissertation, in the design procedure for inductors 

1 

with air-gapped cores a uniform air-gap was assumed. It is suggested that 
analysis of the magnetic circuit with non-uniform air-gap be attempted and 

I 

advantages and disadvantages of inductors with non-uniform air-gap be in- 
vestigated. Thirdly, in the design procedure for air-gapped inductors, 
fringing flux was assumed to be negligibly small. It is suggested that 
design of air-gapped inductors be attempted for the case when fringing flux 
is a significant portion of total flux in the core. 
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Appendix A 

DERIVATION OF EXPRESSIONS FOR THE LC FILTER OPERATING 
IN THE CONTINUOUS-CONDUCTION MODE, CHAPTER II 


The principle of superposition is used to compute steady-state voltages, 
and currents in an LC-fiTter operating in the continuous-conduction mode. The ' 
equivalent circuit of the filter operating in the continuous-conduction mode i 
is shown in Figure 2.2. First, reducing the sinusoidal component of the voltage 

* j 

source v| to zero, the dc currents and voltages in the filter due to the i 

* ! 

I 

application of a dc voltage source of value 2Vj|^/-rr are given as below: i 

2v. i 


V. = — ^ (A2) 

Q IT 

The response of the filter to the second harmonic component of the 
voltage source v| is computed by the sinusoidal steady-state method. The 
filter impedance Z(j2u)) as seen by the sinusoidal source of frequency 2oj 
in radians/sec is given by 


( 

i 


Z(j2o3) 


j2(i)L + 


j2ci)C 

j2wC 


(A3) 


1 

I 
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Eq. (A3) is simplified to 


Z{j2(,) 


(Rj_-4w LCR|^) + j2wL 

l+j2o)CR, 


(A4) 


The magnitude |Z(j2u))l and phase ijj of filter- impedance Z(j2w) may 
be obtained from (A4) in terms of oo.., k, and R, 


|Z(j2(,) 


R|^ ^ 1 -8(ij^+l 6(ijj;^+4K2 




(S/k^) 


ij; = arc tan [ 2 k /(1 -4ai^)] - arc tan (Zw^/k) 


(A5) 


(A6) 


i From Figure 2.2, the value of ij, due to second harmonic component 

f 

is given by 

S 4V, 




Im cos2wt 
3 tt Z(j2(o) 


(A7) 


Substituting |Z(j2(u)j and ip from (A5) and (A6) for Z(j2w) in (A7) leads to 

, _■ 4 ViJl + (4.5/.^) 

1 V ” 


COs( 2 wt-\|;) 


(A8) 


'L 


3ttRi .^1 -8o3m+1 6o)m+4k^ 


N 


Defining 7i(j2to) as the parallel combination of R^ and l/j2wC 


Zi(j2o))l = 


•^l+{4mj^/ic2) 


4> = -arc tan (2 oj|/k) 


(A9) 


(AlO) 
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Thus, the output voltage Vq is given by 


Vq = ij^ Zi(j2co) 


(All) 


Substituting for i^ and Zi(j2o)) from (A8), (A9), and (AlO) in (All) 
and simplifying leads to 

4V, 




Im 


1 -8(ii|+l 6o)J|+4ic^ 


COS(2ojt-ij;+<i)) 


(A12) 


The total inductor current i^^ is obtained by summing the dc and ac 
components given by (Al) and (A8), respectively 


tiR, 


cos{2(jjt-ip) 


(A13) 


iTrRj^ y 1 -8oj^+l 6w|;j+4ic^ 


Using the same approach, the total output voltage Vq is obtained by summing 
components given by (A2) and (A12), 


2V 


v„ = 


Im 


4V 


Im 


0 ir 


3it -^1 -8(jj|+l 6 (u|^+4k^ 


cos ( 2a)t-iii+(j) ) 


(A14) 


From the alternating component in (A14), the peak-to-peak ripple 


i voltage V^p is 


8V 


Im 


'Op 


3-n-^ 1 -8to^+l 6o)^t4K 


(A15) 
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and the average output voltage from (A14) is 


from (A13), the peak and rms values of inductor current, ly and ly , are 


I = 1 + 

^Xm 7 tRj_ 2 


'Xe itR, 


1+(4u!J;/ k^) 


3 y 1 -8 (jo^+1 6o)J^+4k^ 


(1+4(jj|^/k^) 

18(1-8^2+i6io;5+4k2) 


In the case of a center-tapped transformer, current in the two 
secondary windings flows during alternate half cycles and is equal to i,^. 
Thus, the rms value of current in each secondary winding is 


^/’'i 


(A19) 


(A20) 


f 


i Appendix B 

I 

I DERIVATION OF EXPRESSIONS FOR THE LC FILTER OPERATING 

1 

I IN THE DISCONTINUOUS-CONDUCTION MODE, CHAPTER II 

! 

j In this Appendix, the expression for Vq and i^ given by (2.11) and 

I (2.12), respectively, are derived. The derivations of the expressions for 

! . Vn and iy^when the diodes are conducting in the discontinuous-conduction 
[ U a"* 

I 

‘ mode are given in Section B.l. The expressions for the case when the diodes 

I 

\ stop conducting are derived in Section B.2. 


B.l Diodes Conducting 

The equivalent circuit when a diode is conducting is shown in Figure 
2.3A. The state equations for the circuit are 

di„, 

*- “dt ’I 


= Vt - V 


01 


(Bl) 


dv 


01 


dt 


= 1 


'01 


A I 


(B2) 


Differentiating (B2) and multiplying by L leads to 


d^v 


LC 


01 


di 


= L 


XI 


dt‘ 


dt 


L ^^01 
R, dt 


(B3) 


(135) 


136 ' 


Substituting for Ldi^^-j in (B3) from (Bl) and rearranging 


dt 


dt2 \ 0^ I 


(B4) 


Let e = o)t 


dv„T dv^T dvm 

01 ^ 01 de _ 01 


dt de dt de 


(B5) 


^^^01 ^ d ) 2 ^!^ 


dt2 de' 


(B6) 


Substituting for dv^-j and d^Vgi in (B4) from (B5) and (B6) leads to 


dt 


dt2 


'*"''01 ^ U.L '*''01 


„2LC , V,^s1n(9-e„) 


d02 


Rl 


(B7) 


Writing this. last equation in terms of and k 


(jO 


2 l!ol , ‘‘’'Ol 


d02 




(B8) 


B.1 .1 Solution for v 


01 


The particular solution of Vg-j due to the forcing function 
Vii^sin(e-eo) is obtained first. .Assuming the particular solution for v 


01 


to be of the form 


I 
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V = A sine + B cose 

p V V 


(B9) 


and substituting the value of Vp and its derivative from (B9) in (B8) and 
equating the co-efficients for sine and cose from both sides of the resulting 
expression leads to the following set of simultaneous equations for unknowns 
A and B 

V V 






(BIO) 


(BID 


Simultaneous solution of (BIO) and (Bll) yields the following values for 


\ and By 




Vim[coseo(l-w^) -1C sined 


(l-to|)^ + 


(B12) 


B = 

V 


Vim[-sineo(l-w2) -k cosej 


( l - a )|)2 + 


(B13) 


■ The homogeneous part of solution is obtained next. The homogeneous part of 
j the equation (B8) is 


2 ''"''01 . . '*''01 


“N 


de2 


+ K 


V^ = 0 


de '01 


(B14) 


Let the solution be of the form 


"h = 


A0 


(B15) 
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Substituting Vj^ in (B14), the corresponding characteristic equation is 

+ kX + 1 = 0 (816) 

For the discontinuous-conduction case 

4a>§ > <2 (B17) 

Thus, the roots of (B16) are given below 

^1 2 ~ ^ (818) 

where 

6 = -k/2o,2 (B19) 

a = (V ~ ^^)/(2w^) (B20) 

The homogeneous part of the solution for can be written as 

= exp(B0){C^sina0 + D^cosa0} (821) 

The complete solution for Vg^ is obtained by adding the particular 
solution V and the homogeneous solution v, given by-(B9) and (B21), respec- 

r * ' 

tively 

Voi(0) = AySine+B^cos0+exp(60){CySinae+D^cosa0} (B22) 
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From (B21), 


i = c 01 ^ 01 

XI dt R|_ 

(B23) 

or 



(B24) 

Substituting the expressions for Vq-| , and its derivative 
(B24) leads to 

from (B22) into 

\i{e)=F“ A.sine+B.cos 0 -i-exp( 3 e){C-sina 0 +D.cosa 6 } (B25) 

A 1 K|^ 1 I \ 1 

1 where 


(Vi'/'' 

(B26) 

= Bv 

(B27) 

c. = + («|(ec^-.D^))/K 

(B28) 

0, = (“i(SDv-^C^))/ic 

(B29.) 


I B.2 Diodes Not Conducting 

I The equivalent circuit when diodes are not conducting is shown in 

Figure 2.3B. The state equations are 

1X2 = 0 


I 

I 
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f 


(B30) can be written as 


loC 


dv 


02 


'02 


de 


= 0 


Let the solution for (B32) be 


'02 


a e 


A(0-ei) 


Substituting for Vq 2 Tn (B32) from (B33), the root of the 
equation is found to be 


1 _ K 

“CR[_ 


Thus 


Vq2(q) = a [exp{-(e-ei)ic/w2}]. 

at e = Oj, 


From (B35) and (B36) 


'^02 " ^ 02 ^®!^ 


a = Vq 2 { 0 i) 


Thus, 


Vo2(e) = VQ2(ei)[exp{-(e-9i)(c/w2}]^ 


(B32) 


(B33) 

characteristic 


(B34) 


(B35) 


(B36) 


(B37) 


(B38) 


i 


Appendix C 

DERIVATION OF EXPRESSIONS FOR THE LC FILTER WITH NON-ZERO 
SOURCE IMPEDANCE OPERATING IN THE CONTINUOUS-CONDUCTION 

MODE, CHAPTER III 


t 

I 

I 

i 

t 


In this appendix, the expressions for voltages and currents in the 
filter circuit with source impedance included, are derived for the case when 
the filter operates in the continuous-conduction mode. 

From the equivalent circuit in Figure 3.2A and following the method 
in Aopendix A, the dc output voltage Vg and inductor current I^ are given 
by 


Iw = 


ttR|^(p+1 ) 


2V 


V. = 


Im 


0 tt(p+1 ) 


(Cl) 

(C2) 


For computing response to the second harmonic, the filter impedance 

4 

Z(j2aj) seen by the harmonic is 


R. 


Z(j2co) = Rj + j2o,Lgq + 


(C3) 


R,+ - 


L j2wC 


(HI) 
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From (C3), the magnitude [Z(j2a))] and phase ij) of the impedance 
Z{j2co) in terms of 01 ^, sc, p, and is 


|Z(j2o))! = R[_ 


( p+1 -4 (jo^ ) 2 + ( 2sc+( 2pto^ ) /k ) 2 
1 + (4us|5/k2) 


= arc tan [(2 k+(2p(jo^P/k)/(p+1-4o)|)] 
[C2to|)/K ] 


- arc tan 


Following the method in Appendix A, i^ and Vq are given by 


4Vim 


3-rrRj^ ^ (p+1 -4oJ^) ^+( 2 k+{ 2pa)2 )7^)~' 


(C4) 


(C5) 


cos(2wt-!|i) (C6) 


4V 


''o - " 


Im 


3tt (-p+l-4u3^)^H-(2)c+2pa)^/K)^ 


cos(2wt-ijs+<|)) 


(C7) 


Adding the two responses, the total inductor current and output voltage is 
obtained as 


4V 


Iv = 


Im 


X irRj^ I 2 (p+1 ) 




^ l+(4a)!^)/i< 


3 (p+1 -4tC^)2+(2K + 2p(D^j/K)^ 


-cos(2wt-i|j) 


(C8) 
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4V 


^0 II 


Im 


2(p+l ) 


cos(2wt-i|*+4>) 


3 ^ (p+l-4a)^)2+(2ic+2pa3^(/vc)^ 


(C9) 


From (C8), the dc output voltage Vq and peak-to-peak ripple V^p are given by 

2V, 


V = 

0 n(p+l ) 


8V 


Im 


'Op 


3ir ( p+1 ~4w^ ) 2+T2 ic+( 2poj^ )/ k )^ 


(CIO) 


(Cll) 


From (C7), the rms and peak values of inductor current and are ob- 
tained from the expressions 


4V 


I 


Im 


Xe irR, 


(1+(4o)|^P/k:^) 


4(p+l 18( (p+l-4(jj^p2+(2|c+(2ptd|)/K)2) 


(C12) 


4V 


I 


Ira 


Xm irR, 


2lFlT 


4 


(l+(4a}i5)/ic^) 


3 ^ (p+1 -4o)^)2+(2K:+(2pa)^)/K)2 


(C13) 


Appendix D 


DERIVATION OF EXPRESSIONS FOR THE LC FILTER WITH NON-ZERO 
SOURCE IMPEDANCE OPERATING IN THE 
DISCONTINUOUS-CONDUCTION MODE, 

CHAPTER III 


In this appendix, the derivations of the expressions for Vg and 
ij^ when the diodes are conducting in the discontinuous-conduction mode are 
presented. The derivations for the case when the diodes are not conducting 
is the same as given in Section B.2 of Appendix B, and will not be repeated 
here. 

For the equivalent circuit shown in Figure 3.2A, the differential 
] equation for Vq-j can be written following the procedure in Appendix B, 

I Section B.l 


(0 


2LC 


d^v 


01 , ,(oL 


‘‘''01 . ,'*1 


de^ 


. (-.caRjC)-^.(^+l)Voi=Vj^sin(e-eo) (Dl) 


I Writing (Dl) in terms of o)|^, k, and p 


CO 


2 fV 


dv 


01 


de^ 


'r de 


+ (p+l)vg^ = Vj^sin(e-9o) 


(D2) 
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where 


= K + (po)^)/k 


(D3) 


Solution for : 


As described in Appendix B, Section B.l, the particular solution 
for Vqi can be obtained as 


V = A sine + B cose 
p vr vr 


(D4) 


where 


A 


[cose 0 ( P+1 -0)^) - K^sineo] 


vr 


(p+1~u2)2 + J.2 


(D5) 


_ [■ sineo(p+l-wg} - K^coseo] 


vr 


{p+l-0)2)2 + ic^2 


(D6) 


For the homogeneous part of the solution, the roots of the char- 
acteristic equation of (D2) are for the case 4 u) 2 (p+i) > ^2 


'1,2 " ^ 


(D7) 


where 


3^ - -k^/(2oj^P 


(D8) 


cty, - { V “ ic2)/(2ai^) 


(D9) 
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As in Appendix B-2, the solution for Vq^ and i^^^ is found to be 
Voi(e) = By^cose + exp(B^e) [^C^^sina^e + D^^cosa^ej (DIO) 


ixi(e) = 


= A^.^sine + B^.^cose + exp(B^e){C^^sina^e + D^.^cosa^e 


(Dll) 


are 


For the case > 4aj^(p+l), the roots of the characteristic equation 


ai - B^+a' 


(D12) 


“2 = 


(D13) 


where 


a' = 




^r " 


2-1 


The solution for Vq-j and ij^-| is obta-ined as 


(D14) 


where 


Voi(e) = A^v^sine + B^^cose + c;,^exp(ai8) + D;,^exp(a20) 


vr 


vr 


vr 


vr 


(D15) 


ixi(0) = j^A.^sine + B^.^cose + C^lj„exp(ai0) + D^l^exp(a2e)j (D16) 


^ir ^vr 


"N^r 


to^ A 

B, = B +-M-J01 
ir vr K 


(D17) 

(D18) 
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(D19) 


D^r " D^^(l+{o)2a2)/>c2) 

(D20) 

Where 

Dvr> and are given as follows: 


C = 
vr 


{D21) 

O 

< 

II 


(D22) 

c = 

vr 

'*' ^ \r3/(a2"“-|) 

(D23) 

D' = 
vr 


(D24) 


1 

I 

i 

1 

\ 


Appendix E 

DERIVATION OF L. I„^. V^^, f, Aj = 0, CHAPTER IV 

Referring to Figure 4.2, it can be seen that the negative air-gap 
line passes through the point with a slope equal to "Vg^/^g 

and intersects the normal magnetization curve at the point 
negative air-gap line equation can be written: 

Bg - >'o"Sc - 

where N is the number of winding turns and is the permeability of free 

space. Substituting and for H and and rearranging (El) leads 

to an expression for M. 


[ 

i 




Inductance L for the winding, is given by the expression 

y, 

I — ! A ^ ^ \ 

^ ~ + y/,(V^m) 

A g m m 


(E2) 


(E3) 


where y^ is the relative incremental permeability for the core material at 
a given operating point [20], Substituting (E2) into (E3) and multiplying 
the result by S-^^Cy + yields 


(148) 


I 
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which can be rearranged and simplified to 

^^Vo^mQ^mQ^ " ^O^Dc'-^^VmQ^^g 


Equation (E5) can be written in terms of £ ,y and constants as 

g A 


"g + K2 - Ks/y^ = 0 


(E 6 ) 


where K^, K 2 , and K 3 are defined in (4.4), (4.5), and (4.6) respectively. 

Roots of £g for the quadratic equation (E 6 ) in £^ are given by (E7) in which 

only the positive sign for the square root term can lead to a meaningful 

physical value for £ and therefore is shown in the expression. This is 

y 

shown in Appendix 6 . 


£g = 0.5[-Ki + Vk? - 4[K2 - (Kg/y^)]] 


(E7) 


To determine the remaining unknown quantity y^ in (E7), the peak ac 
I flux density must be found. is related to the rms .ac voltage 

across the winding for the assumed sinusoidal fundamental frequency component 
as follows: 




(E 8 ) 


Rearranging (E 8 ) to obtain an expression for N and then squaring leads to 


= K^/B2 


tnAC 


(E9) 
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j where K 4 is given by (4.7). 
i 

Substituting the expression for in (E9) and the expression for 
from (E7) into the equation for L given by (E3) leads to the relationship 
(ElO) which is the desired form Fi(*) = 0. 



j where K 5 and K 0 are defined by (4.8) and (4.9). 


(ElO) 


I 

i 

1 


I 

\ 


i 

i 

1 



Appendix F 

DERIVATION OF THE SCREENING RULE FOR A PHYSICALLY 
REALIZABLE AIR GAP, CHAPTER IV 


The screening process included in Flow Chart 3, Figure 4.7, is based 
on the consideration that, for the air gap to be physically realizable, the 
expression for given by (E7) In Appendix E and in the text as (4.11) must 
yield a real positive quantity. Two sets of conditions will be considered 
by examining the value of the expression for given by (4.4) and repeated 
here as (FI ) 


Ki = mJZZ a H nB n - InrL)/A 
^ 0 m m mQ mQ DC m m 


mQ 


(FI) 


The parameter can also be expressed in terms of K 2 and K 3 , which are de- 
fined by (4.5) and (4.6), and is given by (F2) 


Ki = 2^ - (Kj/i^) 


(F2) 


Case I : Kj 4 0 

Using the definition of in (FI), the inequality 4 0 is equiva- 
lent to the inequality (F3) 




m"m = “mQ^mQ 


(F3) 


( 151 ) 
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and using (F2), it is equivalent to the inequality 

(K3/*„)i2V^ (M) 

The specification-dependent parameter Kg used in the screening pro- 
cess and calculated in Step [I] of Flow Chart 3 is defined as • 


1^8 4 


(F5) 


Thus, for K^ to be equal to or less than zero, (F6) must be satisfied 


Kg > H„nB, 


mQ>Q 


(F6) 


Under this condition, i will be a real number if the discriminant in (E7) 

9 

is equal to or greater than zero or if 


q + (4K3/y^) > 4K^ 


(F7) 


Using the value of K^ from (F2)., the inequality in (F7) reduces to 


(Ka/l,^) + (4yp^) > 4-yfKl (F8) 

Since (F4) holds for the case under consideration, the smallest value for 
I the left side of (F8) occurs with the equality sign in (F4), From (F8), for 
iig always to be a real number with any combination of 5,^, ^ 2 , and Kg, 

then 




(F9) 
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Using the value for K 2 from (4.5) in this expression leads to the inequality 
(FIO) 

^W'o'-Aq " ' <Fio) 

For essentially all quiescent or bias points and all magnetic materials (FIO) 

is satisfied. Therefore, when the inequality expressed by (F 6 ) and shown in 

I Step [J] in Flow Chart 3 is satisfied, the selected core and quiescent point 

I are satisfactory for continuation of the design procedure. 

1 

i 

I Case II: Kj > 0 

For this condition, which corresponds to Ks < order for 

! to be a real positive number, 

I 

I > Kj (Fll) 

> 

j 

I Using the expressions for K 2 and K 3 from (4.5) and (4.6) and rearranging 

j 

, terms leads to inequality (F12) which must be satisfied for a core to be 

I 

! considered further. 

I 

I 

I 

i When the inequality expressed in (F12) and shown in Step [K] in Flow Chart 
3 is satisfied, the selected core and quiescent point are satisfactory and 
the design procedure may be continued. 

If (F12) is not satisfied, a final test is made which compares Ks 
with the value of The parameter Wo(H^P^g)„,.„/2 1s the 


I 
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I minimum value that the product '^^nge of 

! ' 

the magnetic core material data available. If inequality (F13), shown as 
' Step [L] in Flow Chart 3, is satisfied, 

i '^8 <f’ 

Then choosing a new value of H^q in Step [N], which will satisfy the in- 
i equality (F12) in Step [K], will allow the continuation of the design pro- 
’ cedure. If (F13) is not satisfied, a larger core must be selected and the 
screening procedure re-entered. 


I 

! 


i 

I 


I 

I 


Appendix G 

DISCUSSION ON THE QUADRATIC ROOTS FOR i , CHAPTER IV 

y 

As mentioned in Appendix E, only the positive sign in front of the 

square root term is considered in the expression for 5. given by (E7). Vlith 

y 

the negative sign, the expression for becomes 

y 


= 0.5[-Ki - - 4[K2 - (Ka/ii )] ] 


(Gl) 


Case I : Ki 4 0 

Using the definition of Ki in (FI), the inequality Ki 4 0 is equiva- 
lent to the inequality (G2) 


I2 L > 2il A H^nB^n 

DC = m m mQ mQ 


(G2) 


Inequality (FIO) given in Appendix F, which is easily satisfied for essen- 
tially all magnetic materials and bias points, is rearranged as follows 




(G3) 


I Substituting for 2 in (G 2 ) from (G3) leads to inequality (G4) 


iLl > 


'DC' 


0 A m^m“mQ 


(G4) 
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Let 


Kll= 1<2 - (Kg/W^) 


(G5) 


From (4.5) and (4.6) 




*> 0*111 


’ Vmq»A 


(S6) 


From (G4) and (G5) 


Ki3<0 ‘ (G7) 

Expression for in terms of and K^^is 

= 0.5[Ki - 4K,J (G8) 

From (G7) and (G8) it is seen that gap is a negative quantity which is 
not physically realizable, 


Case II - Kj^ > 0 

From (Gl), is either a negative real number or a complex number 
neither of which is physically realizable. 



APPENDIX H 


ALGORITHMS FOR GENERATING DESIGN CURVES FOR AIR- 
GAPPED INDUCTORS, CHAPTER V 


The algorithms used in generating design curves are explained in 

this appendix by reference to the flow chart presented in Figure H.l. In 

the flow chart, starting with an assigned value of air gap ratio smd 

an initial quiescent point the normal magnetization curve, a 

2 

value for design parameter LIqj>/v is computed for an initial value of 

✓s 

®mAC min relationship (5.14). The procedure continues by iteratively 

2 

selecting additional quiescent points until values for LIpj,/v corresponding 
to for all the quiescent points are obtained. The value of 

~ %Q opt quiescent point which corresponds to the 

2 2 
maximum value of LI^^/v, is used to generate design curves of Llg^^/v versus 


B^Ac various values of air-gap ratios using the corresponding H 
for H^q in (5.14). 


mQ,opt 


H.l Flow Chart 

In the procedure illustrated in the flow chart in Figure H.l, a 
search is made for the optimum quiescent point H^q for a specified 
air-gap ratio which yields the maximum value of Llg^^/v. The value for the 

(157) 
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corresponding is used for the dc operating point 

(5.14) to generate design curves of LI^^/v versus shown in Figures 

5.5, 5.6, and 5.7. The operation begins with the assignment of minimum and 
maximum values for the air-gap ratio and the choosing of the magnetic mater- 
ial- The data base for the designated core material includes sample points 
taken from the normal magnetization curve and from a family of relative 
incremental permeability curves for the chosen material, as described in 
Flow Chart 2, Figure 4.. 5, in Chapter IV. 

Starting with the first sample point [Hj^Cl), the normal 

magnetization curve, and using an interpolating polynomial method to ap- 
proximate the values of and for as described in Chapter 

IV, the value for LI^^/v is calculated from (5.14) and stored. This quies- 
cent value is then incremented by a prescribed amount and the above 
procedure is repeated until the specified maximum value for the quiescent 

point is exceeded. The value of H^^g = H^g corresponding to the 

2 

maximum value of LIqj,/v is used now in generating design curve for the 
specified air-gap ratio. 

Using the method of an interpolating polynomial to approximate val- 
ues of B|^g and corresponding to H^g and various values of 
Lljj^/v-is computed from (5.14). A curve of LIqq/v versus B^^^ is obtained 
for the specified air-gap ratio. The value for corresponding 

to H^g ^p^ for the specified air-gap ratio is calculated from (5.15). The 
air-gap ratio is then incremented by a prescribed amount and the 

above procedure is repeated. A family of design curves for grain-oriented 
3-percent silicon steel, obtained in this manner are shown in Figure 5.3. 
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GLOSSARY OF SYMBOLS 
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GLOSSARY OF SYMBOLS 


^0 


A. 




A 


ir 


A„ 


m 


A, 


m,tmn 


A. 

A 


vr 


wn 


wr 

bo to b. 


B 

B. 




ir 


m 


max 


numerical coefficients used in approximate relationships 

given in Table 3.2 and in Table 3.3, numeric 

2 

cross-sectional area of air gap, m 

2 

cross-sectional area of magnetic core material, m 

minimum value of cross-sectional area of magnetic core as 

2 

determined from minimum air-gap volume requirement, m 

[Vijjj(coseo(l-oj^)-Ksin0o)]/[{l-w^)^ + ic^]>V 

[Vj^(cos6Q(p+l-aj^)-K^sin0Q)]/[{p+l-()j^)^ + ic^],V 

2 

area of magnetic core window, m 

2 

cross-sectional area of wire, including insulation, m 

numerical coefficients used in approximate relationships 

given in Table 3.2 and in Table 3.3, numeric 

flux density in a magnetic circuit, T 

flux density in air gap, T 

(A^w^/k) + V 

^ ®vr- '' 

flux density in magnetic material, T; corresponds to 
values on the normal magnetization curve 
maximum allowable core flux density, T 
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166 


"mAC 


B, 


mAC,max 


B, 


mAC,min 


peak value of the fundamental component of the ac flux 
density in magnetic material, T 
maximum value of B^ 

termined from incremental permeability data, T 


for a particular value of as de- 


minimum value of B^^^ for a particular value of as de- 
termined from incremental permeability data, T 


"mDC 


"mQ 


B, 


mQjOpt 


vr 

c^ to c^ 




ir 

CU 


vr 

C 

vr 


d-| to d^ 




“ir 


average or direct component of B^^, T 

selected quiescent or bias-point value of B^, T 

value of Bj^q corresponding to optimum (minimum) turns for 

a specified value of inductance, T 

residual core flux density, T 

-[Vi|^{sin0o(l-(u^) + KCOseQ)]/[(l-u^)^ + k^], V 

-[Vi|^(sineo(p+l-w^) - K^cos6g)]/[(p+l-w^)^ + V 

numerical coefficients used in approximate relationships 

given In Table 3.2 and in Table 3.3, n umeri c 

capacitance value of the output capacitor, F 

(a)^(BCy - aDy)/K) + C^, V' 

(co^(B C - a D )/k) + C , V 
\ N r vr r vr*^ ‘ vr 

Cy^[((u^a.j)/K) + 1], V 

symbol used to identify current stepup power stage 
[(-KVj^sin0q/w^) - Ay - 3(Vjj^sin0o-By)]/a, V 
[(-KcVj^sineQ/to^) - Ay^ - B,(Vj^sin0Q-By^)]/a^, V 

^“2^'^Im^^'^®0“^vr^ + Ay^]/(a2-a.j ) , V 

numerical coefficients used in approximate relationships 

given in Table 3.2 and in Table 3.3, numeric 

w^{BDv+aCy)/K + By, V 

o 

cOm(B D +CX C )/k + D , V 
r vr r vr' vr 
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ir 

= 

Dy^[((w^a^)/K) + 1], V 

1 

V 

= 


1 

vr 

= 


i' 

vr 

= 



frequency of the fundamental component of voltage across 
the reactor, Hz 

symbol used to identify constant frequency controller 

winding factor, numeric 

maximum permissible winding factor, numeric 

magnetizing force. A/m 

magnetizing force in air gap, A/m 

magnetizing force in magnetic material. A/m; corresponds to 

values on normal magnetization curve 

average or direct component of A/m 

selected quiescent or bias-point value of A/m 

value of corresponding to optimum (minimum) turns for 

a specified value of inductance. A/m 

maximum of the average or direct component of inductor 

current, A 

effective or rms value of input current, A 
instantaneous value of output current, A 
instantaneous value of output current when diodes are 
conducting, A 

instantaneous value of output current when diodes stop 

conducting, A 

average output current, A 



168 


^Pe 

^Pm 

''Sl ?^S2 

'Se 

^Sm 

^'x 

'X1 

^X2 

^Xe 

^Xe,max 

^Xm 

K4 

^5 

^6 

^7 

^8 


instantaneous value of current in the primary winding of 
transformer, A 

effective or rms value of current in the primary winding 

of transformer, A 

peak value of primary current, A 

instantaneous value of currents in one of the secondary 
winding of the transformer, A 

effective or rms value of current in each secondary wind- 
ing of the transformer, A 
peak value of secondary current, A 
instantaneous value of inductor current, A 
instantaneous value of inductor current when diodes are 
conducting, A 

instantaneous value of inductor current when diodes stop 
conducting, A 

effective or rms value of inductor current, A 
maximum value of effective, inductor current, A 
peak value of inductor current, A 

^O^m^OC^'^VmQ’ ^ 

(VAc/V^TrfAj2, 

21 > m 
m 

L/2yQA^, m"^ 

-K^ - 4[K2 - , m 

InpL/2A A , HA^/m^ 

DC mm 
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numeric 


(V2Po.«x»''l.™in-''())<B„„-B,)/(Vo^V„-V(,),VU.(AT)- 




'V2Po.max)(''l,m1n-''q»8«,ax-8R)''( 1 ,^uV*W 


772 l^'‘o‘'AWniQ ■ ^ 

VmQ“A 

length of air gap, m 

minimum value of air-gap as determined from minimum air- 
gap volume requirements, m 

length of air gap corresponding to optimum (minimum) turns 
for a specified value of inductance, m 

mean magnetic path length in the magnetic core material, m 


inductance value of inductor, H 

effective inductance associated with input voltage source 
as referred to the secondary side of the transformer, H 
inductance of primary winding, H 
inductance of secondary winding, H 

harmonic number in the fourier series expansion of input 
voltage, numeric 

number of turns on inductor winding, numeric 
minimum number of turns on inductor winding for a specified 
value of inductance, using a given core material, geometry 
and air gap-, numeric 

number of turns on the primary winding, numeric 
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P.F. 

p 

0,max 

Q 

\ 

t 

^off 

^on 

T 

TF 

TN 

"y 

UD 

2UD 

''ac 

''o 

''l 

''le 


- number of turns on each secondary winding, numeric 
= power factor of the rectifier-filter circuit, numeric 
= maximum value of average output power, W 
= quality factor of the filter at resonant frequency, numeric 
= effective resistance associated with input source as re- 
ferred to the secondary side of the transformer, 

= output load resistance, 

= time, s 

= transistor cut-off interval, s 

= transistor conduction interval, s 

= period of one conversion cycle, s 

= symbol used to identify constant off-time controller 
= ■ symbol used to identify constant on-time controller 
= option constraint quantity, y = A, B, C, D, E, H, or J; 
units depend on option 

= symbol used to identify single winding voltage stepup/ 
current stepup power stage 

= symbol used to identify two winding voltage stepup/current 
stepup power stage 

= minimum of the rms or effective value of the fundamental 
alternating component of voltage across the inductor, V 
= diode forward voltage drop, V 

= converter DC input voltage, V 

«• 

= effective or rms value of ac input voltage to the rectifier- 
filter circuit, V 

= peak value of the sinusoidal input voltage, V 
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''O 

''oi 

'^02 

Vq 

^Oe 

'op 

% 

'q 

vu 

|Z(j2w)| 

lZ^(j2w)i 


a' 

“l 

“2 

“r 


maximum converter input voltage, V 
minimum converter input voltage, V 
instantaneous output voltage due to homogeneous part of 
the differential equation, V 
instantaneous value of output voltage, V 
instantaneous value of output voltage when diodes are con- 
ducting, V 

instantaneous value of output voltage when diodes stop 
conducting, V 
dc output voltage, V 

effective or rms value of output voltage, V 

peak-to-peak output ripple voltage, V 

instantaneous output voltage due to forcing function of 

differential equation only, V 

transistor saturation voltage drop, V 

symbol used to identify voltage stepup power stage 

magnitude of the filter impedance seen from the source at 

the second harmonic of the sinusoidal ac supply, fi 

magnitude of the parallel combination of filter load and 

capacitance at the second harmonic of the sinusoidal ac 

supply, 

,/4oj^-ic^/(2oj^) , numeric 

,/{icJ-4u^(p+l)/{2w^), numeric 

3^ + a', numeric 

3^ - a' , numeric 

I ? 2 

V4o)^(p+l) - numeric 
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g 

Y 

6 


^1 

e 

^max 

^min 

C 

6 


X 

^DC 


2 

-k/(2o)|^), numeric 
2 

-k^/(2w^), numeric 

Ng/Np = turns ratio, numeric 

turns ratio N^/Np for option y where y = A,B,C,D,E,H, or 0, 
numeric 

computed quantity used in calculating minimum air-gap 
volume (Reference 12), m^/H 

error term computed in solution for discontinuous mode, 
numeric 

error term computed in solution for and numeric 

/s 

value of e corresponding to and numeric 

value of e corresponding to and numeric 

computed quantity used in calculating 6, N, and Np (Refer- 
ence 12), m"^/H 
wt, radians 

delay angle between beginning of sine pulse and resumption 
of inductor current in the discontinuous mode, radians 
end of the conduction interval in the discontinuous mode, 
radi ans 

toL/R|^, numeric 
2 

(c+(pWj^/k), numeric 

critical value of the parameter k, numeric 
coefficient of the exponent in the solution for the dis- 
continuous mode, numeric 

relative normal permeability of the magnetic material, 
numeric 



173 


V,eff 

^eff 

^0 

^A,eff 

^A,max 

^AjFiiin 

^AQ 


V 

V . 

tmn 

? 

P 

4> 


effective relative normal permeability of the air-gapped 
structure, numeric 

effective permeability of the air-gapped structure, 
numeric 

permeability of free space, 4 tt x 10 ~^ H/m 
relative incremental permeability, numeric 
effective relative incremental permeability of the air- 
gapped structure, numeric 

maximum relative incremental permeability corresponding 

to B .p numeric 
mAC,max 

minimum relative incremental permeability corresponding to 

A 

numeric 

mAC,min 

A 

^A^%’^mAC,max^ "" maximum value of on relative incre- 
mental permeability curve corresponding to the selected 
quiescent or bias point Hj^q, numeric 

3 

I = volume of magnetic core, m 
mm ^ 

minimum volume of the powder permalloy or ferrite toroidal 
3 

cores, m 

damping ratio, numeric 

Rj/Rj^ = ratio of source resistance Rj to output load re- 
sistance Rj^, numeric 

phase of the equivalent impedance resulting from parallel 
combination of filter load and capacitance at the second 
harmonic of the sinusoidal ac supply, fi 
phase of the filter impedance seen from the source at the 
second harmonic of the sinusoidal ac supply, radians 
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ca 


to 


0)k 


Ae 


AH, 


m 


g m 


AW, 


m.max 


27rf, angular frequency of the input sinusoidal source, 
radians/sec 

1/ LC, the natural angular frequency of the filter with 
zero damping, radians/sec 

u)/caQ» normal i zed angular frequency of the sinusoidal input 
source, numeric 

range factor for comparison of error term e, numeric 
incremental change in used in optimization algorithm to 
determine minimum turns for a specified value of induct- 
ance, A/m 

incremental change in (^g/^^) in the algorithm for gener- 
ating design curves, numeric 

maximum amount of energy transferred by core over a 
switching cycle, 0 
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